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Supersymmetry searches at the LHC are both highly varied and highly constraining, but the 
vast majority are focused on cases where the hnal-stage visible decays are prompt. Scenarios 
featuring superparticles with detector-scale lifetimes have therefore remained a tantalizing 
possibility for sub-TeV SUSY, since explicit limits are relatively sparse. Nonetheless, the 
extremely low backgrounds of the few existing searches for collider-stable and displaced new 
particles facilitates recastings into powerful long-lived superparticle searches, even for mod¬ 
els for which those searches are highly non-optimized. In this paper, we assess the status of 
such models in the context of baryonic R-parity violation, gauge mediation, and mini-split 
SUSY. We explore a number of common simplihed spectra where hadronic decays can be 
important, employing recasts of LHC searches that utilize different detector systems and 
hnal-state objects. The LSP/NLSP possibilities considered here include generic colored su¬ 
perparticles such as the gluino and light-flavor squarks, as well as the lighter stop and the 
quasi-degenerate Higgsino multiplet motivated by naturalness. We hnd that complementary 
coverage over large swaths of mass and lifetime is achievable by superimposing limits, par¬ 
ticularly from CMS’s tracker-based displaced dijet search and heavy stable charged particle 
searches. Adding in prompt searches, we hnd many cases where a range of sparticle masses 
is now excluded from zero lifetime to inhnite lifetime with no gaps. In other cases, the 
displaced searches furnish the only extant limits at any lifetime. 



I. INTRODUCTION 


In spite of progressively tightening limits 0121 , supersymmetry (SUSY) continues to serve 
as one of the most compelling scenarios for new physics at the Large Hadron Collider (LHC). 
In the absence of traditional signals below about 1 TeV, there has been a growing interest 
in exploring models in which superparticles are light enough for the 8 TeV LHC to produce, 
but somehow manage to evade the existing searches. Indeed, this situation remains a strong 
possibility, in large part because the observable signatures of SUSY depend extremely sen¬ 
sitively on the details of the spectrum and on the decays of the lightest superparticles. In 
particular, a broad class of highly-motivated scenarios that has so far received relatively 
limited dedicated attention in LHC searches includes superparticles with macroscopic decay 
lengths spanning from sub-mm to tens of meters. Such “displaced” particles can occur in 
models with R-parity violation (RPV) [3] or gauge mediated SUSY breaking (GMSB) |1], 
and also in mini-split spectra where all squarks are at roughly the 1000 TeV scale EE] and 
the gluino lifetime becomes extended. 

However, displaced decay signals occupy a subtle place in collider phenomenology. At 
the high end of the above lifetime range, long-lived charged particle searches may 

become appropriate if the displaced particle is charged or has a nontrivial chance to form 
charged hadrons, but the efficiency drops exponentially for lower lifetimes. At the low 
end of the lifetime range, any number of prompt searches or searches involving bottom or 
charm hadrons might pick up the signals Hi, but may subject them to unnecessarily large 
backgrounds. Many different targeted strategies have been applied to search the broad range 
lifetimes in between [TT]j21]- Because energetic particle production originating in the bulk 
of the detector is extremely rare, most of these searches beneht from tiny backgrounds, 
often 0(1) event or smaller while maintaining good signal efficiency. The non-observation of 
excesses in such clean searches then begs the question of their implications for more general 
classes of models, where even very low efficiency can still lead to signihcant limits. Several 
recent phenomenological works have investigated the power of these searches, or proposed 
new searches for similarly striking signals [2SH3Zj. 

In this paper, we attempt to develop a more rehned understanding of the status of SUSY 
scenarios with displaced decays, in light of the small but powerful collection of existing 
LHC displaced particle searches. Most of the models that we study have either never been 
searched for at nonzero lifetime, or explicit searches cover only one possible signature out of 
several. In practice, it is typical for each highly specialized displaced decay search to phrase 
its results in terms of only one or a handful of highly specialized new physics models. This is 
understandable, given the vast range of possible interesting models and the computational 
overhead required to fully simulate and interpret them. However, inferring the implications 
of those searches for a different model, and in particular how they might interplay with 
each other in covering the parameter space of that different model, then requires careful 
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recastings. We present here, what we believe to be for the hrst time in the context of our 
chosen SUSY models, a comprehensive set of such recastings for multiple displaced particle 
searches simultaneously. Our results highlight the strengths and weaknesses of the various 
searches, and give a clear indication of what regions may currently be lacking in sensitivity. 
The process of undertaking these recastings also illustrates some of the difficulties and 
ambiguities that can arise when attempting to extrapolate the results of displaced decay 
searches beyond their original target models, especially given the unconventional approaches 
to event reconstruction. A saving feature is that total rates near the boundaries of sensitivity 
are usually very strong functions of both mass and lifetime, such that even 0(1) uncertainties 
in our estimation of experimental acceptance can still lead to only 0(10%) uncertainties in 
model reach. Nonetheless, where possible we point out aspects of the searches that could be 
particularly prone to mismodeling by recasters. We also make several suggestions for how 
some of these searches might be adapted to serve as more powerful probes of SUSY or other 
models beyond their original targets. 

Of course, even restricting ourselves to simple variations on the particle content of the 
Minimal Supersymmetric Standard Model (MSSM), the variety of possible displaced hnal- 
state signatures is extremely rich. To narrow down the possibilities to a manageable level, 
we hrst of all focus on simplihed models where only one type of superparticle is produced, 
and undergoes a single-stage displaced decay back into SM particles and/or the LSP. These 
simplihed models can generally be embedded into a variety of more complete spectra with 
additional production channels, such that our results are both broad in applicability and 
conservative within any given model. (For a more inclusive survey approach, see [35].) 
Within the still rather large set of possible simplihed models, we focus on ones that have 
a sizable fraction of hadronic visible decays, either directly or due to subsequent decays of 
electroweak bosons. Such hadronic signals are nominally the most challenging and the least 
constrained by explicit displaced searches, and in some cases unconstrained or only mildly 
constrained even in prompt decay searches. Signihcantly, some of the simplihed models 
that are most motivated by naturalness [3B1 - IT2] can fall into this category, including direct 
production of the lightest stop eigenstate or of a quasi-degenerate multiplet of Higgsinos. 
All together, the models that we consider include: 

• t —>■ didj via baryonic RPV, including i ^bb [SU] (Figs. 

• g Uidjdk via baryonic RPV (Fig. 

• H ^ Uidjdk (-|-soft) via baryonic RPV (Fig. 

• q^qG in GMSB (Fig. g 

• g^gG in GMSB (Fig. Q 

• G in GMSB (Fig. g 
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• H ^h/ZG (+soft) in GMSB (Fig. § 


g —>■ qqB in mini-split SUSY (Fig. 10) 


Some major options missing from this list are sleptons, electroweak gauginos, simplified 
spectra with leptonic RPV, and mini-split SUSY with gluino decays dominated by heavy 
flavor. As discussed in more detail below, some of these other possibilities are covered 
already by existing searches or recasts, and some we expect to have signihcant overlap with 
the above signals, but some would also be worth a closer look in future work. 

The most powerful displaced decay limits within our selection of models typically come 
from the CMS tracker-based search for displaced dijets HT]. This search often remains 
sensitive to models with cr much larger than the tracker radius, as well as to models with 
decay topologies different from the nominal dijets. For models where the long-lived particle 
is colored, hadronization implies a sizable charged fraction that can also be picked up by 
stable charged particle searches in events where the decay takes place outside of the detector. 
Similarly, these searches maintain some sensitivity for cr much smaller than the 5-10 m 
outer detector radius. The overlap of exclusions between displaced decay searches and 
stable charged particle searches can then be signihcant, sometimes more than three orders 
of magnitude in lifetime. At the low end of the lifetime range, prompt searches also become 
sensitive. While it is not possible for us to precisely map out the lifetime range over which 
these searches remain efficient, conservative guesses again allow for signihcant overlap. This 
complementarity often allows for exclusions that span from prompt lifetimes to inhnity with 
no gaps. 

The main results of this paper consist of a series of exclusion plots over the mass-lifetime 
plane of each displaced particle. Figs, [^through For the colored production models, 
the mass reach in the cr range of O(mm-m) is usually comparable to, and in some cases 
better than, the ~ 1 TeV reach from collider-stable charged particle searches. In particular 
stops, which are expected to have mass less than about 1 TeV in a natural model, have very 
little viable model space surviving in this lifetime range under these decay scenarios. For the 
electroweak Higgsino production, stable charged particle limits do not apply, and prompt 
searches are typically limited in sensitivity, but a large number of displaced searches yield 
powerful limits, especially in the GMSB case. We hnd that for cr ~ 10 cm, masses below 
about 600-800 GeV are excluded, giving serious tension with naturalness at those lifetimes. 
For natural masses near 100 GeV, the excluded lifetime ranges from 0(10 microns) to 
0(10 m) in RPV, and up to 0(100 m) in GMSB, dominated there by GMS’s tracker-based 
displaced dilepton search ng. 

For all models, the region of lifetimes around cr ~ 10 m could in principle beneht 
from searches in the hadronic calorimeters and muon chambers, such as those performed 
by ATLAS [121 12D]- But the existing searches are highly limited in sensitivity by their 
focus on lower-mass models and by requiring very tight reconstruction cuts on both sides 
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TABLE I: A summary of the LHC searches recast in this paper. 


of the event. For the Higgsino models, improvements in this direction might be the only 
option for extending the sensitivity to higher lifetimes, without ultimately appealing to 
more standard-style SUSY searches that assume that assume that both hnal-state Higgsinos 
escape the detector unseen. 

The paper is organized as follows. In the next section, we review the existing LHC collider- 
stable and displaced particle searches that we use in our limit-setting. (This section may 


be bypassed by a reader who is not interested in the details of these analyses.) Section III 


specihes the motivations and features of the simplihed SUSY models under investigation, 
and presents our derived limits. We conclude and present some ideas for future searches in 


Section ]^ An appendix discusses the details and calibrations of our detector simulations 
used for recasting. 


II. THE LHC SEARCHES UNDER CONSIDERATION 

Displaced decay searches at the LHC are currently limited to a handful of specihc new 
physics scenarios.^ Searches that target minimal SUSY include non-pointing photons in 
gauge mediation (assuming a mostly-bino LSP) [II1II2], fhe “disappearing track” signature 
of NLSP charginos in anomaly mediation [T31135] . displaced leptons from neutralino or stop 
decays with leptonic RPV caiisiEi, and late decays of gluino R-hadrons stopped in the 


^ Displaced decay searches have also previously been carried at the Tevatron [TlllHT] and at LEP (STHSU- 
These searches have for the most part either been superseded by the LHC or do not have immediate 
relevance to the SUSY models we consider. We do not attempt to recast any of them. However, we 
practically assume that long-lived particles below 100 GeV should have been highly visible to some of 
these searches. 
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calorimeters in mini-split SUSY [T6l Other searches focus on models such as Hidden 
Valleys [171120112311251126] or light hidden-sector gauge bosons [211 122] • Recently, ATLAS has 
also re-interpreted its prompt gluino limits, accounting for the effect of displacement on the 
signal acceptance [lO], results that we put into broader context here. CMS has re-interpreted 
its stable charged particle searches for a large ensemble of pMSSM models with long but 
hnite lifetimes [0], and we apply a similar strategy to our more focused set of models. 

From this modest but growing list of analyses, we select seven that appear to be of greatest 
relevance for the SUSY models studied in Section the CMS heavy stable charged particle 
search |8|, the CMS displaced dijets search HZI, the CMS displaced dileptons search [T8] . 
the CMS displaced electron and muon search na, the ATLAS muon chamber search [20]. 
the ATLAS low-EM jet search [T9|, and the ATLAS displaced muon plus tracks search na. 
Except for the ATLAS muon chamber, all of these have been performed with the full 8 TeV 
dataset. The following subsections summarize the relevant aspects of each analysis that we 
use for our recasts, as well as commentary on the reliability of these recasts where appro¬ 
priate. Our approximate reproduction of each of these analyses relies on simplihed detector 
simulations. Descriptions of these simulations and their calibration to known experimental 
results is provided in a corresponding set of subsections in the appendix. Table [I] provides a 
compact overview, including the associated references and subsections. 


A. CMS Heavy Stable Charged Particles 

Both ATLAS [7] and CMS [H] have conducted searches for heavy stable charged parti¬ 
cles (HSCP) that traverse the entire detector, and appear as a “heavy muon” with anoma¬ 
lously small velocity or dE/dx in the detector material. We choose to focus on the CMS 
searches [8], though we expect very similar performance from the ATLAS searches. Stable 
squarks and gluinos have already been explicitly considered by both experiments, and will 
simply generalize these results to cases with finite lifetimes. 

The only major subtlety when dealing with meta-stable colored particles is that they are 
only seen bound into R-hadrons isg. The hadronization fractions can be estimated from 
simple models, and are likely fairly accurate for squarks given the extensive theoretical and 
experimental experience with heavy quarks. Hadronization of the color-octet gluino is less 
certain, but we assume here the default behavior in PythiaS [58] • (This results in a charged 
hadronization fraction of approximately 46%.) A more subtle issue is how these R-hadrons 
interact with the detector material, especially the chance that a charged R-hadron will pass 
through the calorimeters without a net charge exchange, and thus manage to trigger in the 
muon system. CMS considers two models: a nominal hadronic cloud interaction and a more 
extreme “charge-stripped” assumption where all R-hadrons emerge from the back of the 
calorimeter in a neutral state. Thankfully, the complicated interplay with the detector has 
been accounted for by CMS, and to extract our own finite-lifetime limits we can concentrate 
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on simpler, geometric considerations. 

We consider two of their search strategies. For the nominal hadronic interaction model, we 
take the tracker plus time-of-flight analysis. For the pessimistic charge-stripped assumption, 
we take the tracker-only analysis. 

The tracker plus time-of-flight analysis relies dominantly on the muon trigger, and 
searches for anomalous track candidates that are matched between the muon chamber and 
inner tracker. This track must be reconstructed with \ri\ < 2.1 and pt > 70 GeV, inverse- 
velocity above 1.225/c (measured using timing information), and a high dE/dx. There are 
also additional requirements on the mass inferred from the momentum and velocity measure¬ 
ments, which are constructed to be highly efficient for signal. In order to recast the cross 
section limit for a given model at hnite lifetime, we form a conservative rescaling factor. 
The numerator is the number of charged R-hadrons that would pass into the analysis given 
the above cuts, excepting the dE/dx cut, which we cannot model but which should also 
be highly efficient for signal. The denominator is the number of charged R-hadrons that 
pass these criteria and that also decay fully outside the detector. We only consider decaying 
R-hadrons where none of the visible decay products re-intercept the detector volume, as this 
may cause additional activity in the muon chambers and could have an unpredictable effect 
on the acceptance.^ 

In order to access the charge-suppressed scenario, the tracker-only analysis exploits a 
subtlety of the trigger. A charged R-hadron may leave a track in the inner tracker, but if 
it leaves no track in the muon chamber and minimal calorimeter activity along its trajectory, 
the particle-flow algorithm used in triggering will assume that the track is spurious and not 
count it toward the calculation. The R-hadron therefore adds to the apparent 
Because each R-hadron either leaves such a “trigger-invisible” track or is neutral to begin 
with, the apparent is the total recoil pt of the heavy particle pair.^ Offline, events 
from this .^T-triggered sample can be analyzed for inner tracks consistent with heavy stable 


^ CMS has also provided a full efficiency map of this analysis [S] , which can be extremely useful in general 
recasts. However, we do not use this map since our physics models are identical to the ones that CMS 
studies, up to the finite lifetime. There could in principle be some interplay between the variation in 
stable particle acceptance and non-decay probability versus kinematics, which we are not simulating, but 
our treatment should be conservative. For example, slow particles would tend to decay earlier and become 
vetoed from the analysis, but slow stable particles (especially with jS < 0.4) are anyway less efficiently 
accepted. Similarly, particles at higher |? 7 | must survive over a longer three-dimensional path length before 
exiting the detector, and again are more likely to be lost due to decay, but high-|ry| is also less efficient 
even for stable particles. Therefore, our naive approach, which effectively assumes a flat acceptance within 
the fiducial region, misses the fact that the particles that survive undecayed also tend to be in kinematic 
regions with higher acceptance. In any case, the turn-off of overall acceptance for this analysis due to 
decays at lower lifetimes is exponential, and we expect this behavior to dominate. 

^ For all SUSY pair production, processed through PythiaS, we damp the ISR, which has been shown to 
better-reproduce matched results [55] . 
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particles. The basic track |? 7 | and px requirements are the same, but there is no velocity cut 
(as no timing information with respect to the muon chambers is available), and the dE/dx 
requirement is tightened. Again, we cannot model the dE/dx cut, so we assign an ad hoc 
velocity ceiling of 0.7c, which puts us on the steep section of the Bethe-Bloch stopping power 
curve [BD]. (Our final results are not very sensitive to the placement of this velocity cut.) We 
again form a rescaling factor for the inhnite-lifetime cross section limits presented by CMS. 
For the numerator, we take the number of events where the recoil px exceeds 150 GeV and 
at least one R-hadron is charged and passes the reconstruction cuts. The denominator is the 
number of events that satisfy these criteria and where both R-hadrons decay fully outside 
the detector (including the non-intercept requirement on the visible daughters, as above). 


B. CMS Displaced Dijets 

The CMS displaced dijet search [T71 I3B] uses a specialized trigger to capture events 
containing a pair of high-p-p jets containing displaced tracks at the level of several hundred 
microns. In the offline analysis, it counts the total number of jet-pairs that appear to be 
consistent with common vertices with a large number of such displaced tracks. For our 
recasts, we focus on the “High-Lxy” analysis, which has 1.14 ± 0.54 expected background 
vertices and one observed, placing an upper limit of 3.7 signal vertices. We have found 
that the High-Lxy works well for all of our models, even ones with short decay lengths 
((Lxy) < 20 cm), and that the choice of High-Lxy versus the very similar Low-Lxy analyses 
has only minor impact on our results. 

The jets used in the analysis have px > 60 GeV and |p| < 2, and the total event Hx 
must exceed 300 GeV at trigger level. (In practice we use a slightly tighter 320 GeV to 
account for the observed turn-on of the trigger with Hx measured offline |6T].) Each pair 
of such jets is inspected for associated tracks with impact parameters larger than 500 /im, 
and this set of displaced tracks is checked for consistency with a common vertex. At least 
one displaced track from each jet is required to fit that vertex, and there are a number of 
additional quality requirements on the vertex itself: total track-mass greater than 4 GeV, 
total track-pT’ greater than 8 GeV, “significant” transverse distance Lxy from the primary 
vertex, and a multivariate likelihood-ratio discriminant cut. The discriminant is formed from 
distributions over the vertex track multiplicity, the fraction of tracks with positive impact 
parameters (based on the sign of the dot product of the track’s px vector and transverse 
displacement vector at the transverse point of closest approach to the beamline), and two 
additional variables based on a special clustering of track crossing points along a line starting 
at the detector’s center and oriented with the dijet px direction. The exact algorithm for 
this clustering is not given by GMS. For our reproduction of the analysis, we create a 
sliding window of full-width 0.15 x Lxy, and adjust it to surround a maximal number of 
crossing points. When multiple window locations would surround different crossing point 
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collections with equal multiplicities, we choose the one with the smallest RMS. For its 
multivariate discriminant, CMS includes both the cluster multiplicity and the RMS relative 
to the vertex Lxy. We determine the event-by-event value of the discriminant by using CMS’s 
own distributions for the four variables. (Our calibration distributions for these variables 
and the multivariate discriminant can be found in Appendix]^) 

The rest of the High-Lxy selection demands at most one prompt track per jet within the 
dijet candidate, less than 9% of the energy of each jet associated to prompt tracks, and 
a multivariate discriminant value greater than 0.8. In the analysis note [56] (which has 
identical results as the more recent preprint CZI). CMS’s new physics limit is phrased in 
terms of the number of dijet candidates that pass all of these requirements over the full 
8 TeV run. 

While this defines the basic search, we point out a few possible subtleties: 

• CMS has only performed full simulation on models with dijet masses up to 350 GeV. 
We assume that there are no major obstacles to probing masses beyond 1 TeV. These 
may experience more tracking confusions due to the greater multiplicity of hits, but 
on the other hand should also be capable of surviving the analysis cuts with a smaller 
fraction of successfully-reconstructed tracks. 

• For all of the physics models that CMS has studied, the displaced particle is neutral 
and leaves no tracker hits before its decay. It is unclear what would happen for charged 
displaced particles, such as R-hadrons, which would leave a signature sometimes called 
an “exploding track.” Presumably the extra hits would lead to additional confusions 
of the tracker algorithms. We simply exclude such cases from the analysis, effectively 
assuming zero efficiency. This is certainly an over-conservative treatment, especially, 
for decays that occur before crossing the first pixel layer. 

• For some of our SUSY models, more than two jets can originate from the same dis¬ 
placed vertex. In principle, in the note [56] CMS considers all possible displaced dijet 
pairs, allowing the same jet to appear multiple times. Such a decay could therefore 
contribute much more than one dijet candidate. However, since the exact procedure 
is not unambiguously described in the CMS note, when a dijet pair passes the basic 
selection cuts (before the High-Lxy selections), we remove its constituent jets from 
further consideration for constructing other dijet pairs. 

• For most of our SUSY models, the displaced dijet candidate will be produced in a decay 
with additional activity, and will therefore by itself not reconstruct the decaying par¬ 
ticle’s momentum vector. Since CMS’s cluster discriminant variables are constructed 
under the assumption that the vertex displacement vector and dijet momentum vector 
are well-aligned, there is a question of whether the discriminant is particularly inef¬ 
ficient for models where this is no longer true. We have found that any such effect 


is quite minor, and that the multivariate discriminant is mainly driven by the vertex 
track multiplicity variable.^ 

• When a decay contains bottom and/or charm quarks, it may generate multiple nearby 
displaced vertices rather than a single displaced vertex. CMS gives some explicit 
indication of how the reconstruction rate differs for heavy flavor decays, and it appears 
that such small secondary displacements do not play a signihcant role, but only for 
models down to 0(cm) lifetimes. The behavior for shorter lifetimes is not specihed. 
Due to the ambiguity, we simply ignore decays that contain heavy flavor and have 
Lxy < 1 cm. 

• It is clarihed in the more recent analysis preprint lizi that no lepton identihcation 
is utilized, and that electrons and muons would function as jet constituents in this 
analysis. However, the earlier note [56], on which we base our analysis, is not explicit 
on this point. We have therefore excluded isolated leptons (as per the dehnition of |18j l 
from jet clustering. This may lead to slightly over-conservative limits when leptonic 
decays are available. 


C. CMS Displaced Dileptons 

The CMS displaced dilepton search [18] is in some ways a simpler version of the displaced 
dijet search above.® The analysis operates on pairs of isolated or , demanding that 
a pair reconstruct to a common highly-displaced vertex. The analysis is both background- 
free and has zero observed events, leading to a 95% limit of approximately 3 signal events. 

Electron-pair candidates in the event must consist of a leading (subleading) electron with 
Pt > 40 GeV (25 GeV). Muon-pair candidates must consist of muons with pt > 26 GeV. 
Each lepton should be in the well-instrumented region of the tracker, at \t]\ < 2.0, and should 
have a “signihcant” transverse displacement of 12cr relative to the resolution. Since here we 
do do not simulate track-by-track resolutions, we simply replace this last criterion with a 
hxed cut of 250 pm. This choice appears to give conservative results for very low lifetimes. 
GMS further requires the leptons to be isolated from other tracks, excluding other identihed 
leptons, using hollow cones of outer radius 0.3 built around the lepton candidates. Since our 
own lepton ID is “perfect” (up to the efficiency factors discussed in the appendix), we use 
solid cones. This again should furnish a conservative approximation, and should roughly 
approximate the ID failures that would occur in reality, e.g. for leptons inside of 6/c-jets. 
Non-lepton tracks with px > 1 GeV count toward the isolation, and should tally to less 
than 10% of the lepton px- For muon-pairs, the two candidate tracks must be separated 


^ The same conclusion was reached in |M| . 

® We do not consider the superceded version of this search at lower luminosity and beam energy |62j . 
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by AR > 0.2, to operate in the region of high dimuon trigger efficiency. For electron-pairs, 
we here additionally require AR >0.1, assuming that the ECAL patterns of closely-spaced 
electrons could start to become difficult to reconstruct. A number of other basic quality cuts 
are applied: the dilepton pair should have invariant mass above 15 GeV to avoid hadronic 
resonances, the azimuthal angle between the dilepton momentum vector and its displaced 
vertex position vector should be less than 7r/2, and the 3D opening angle between a pair of 
candidate muon tracks should have a cosine greater than —0.79 to veto cosmics. 

Given the simpler topology and cuts relative to the displaced dijet search, we should 
naively be more confident in the robustness in our recasting procedure for different models. 
Nonetheless, as discussed in Appendix]^ our modeled acceptances appear to be low for some 
of GMS’s benchmark models, possibly because we err on the conservative side for the lepton 
track-finding efficiencies at large decay radii. The “exploding track” question also persists 
in principle, but does not arise in any of the models that we consider. (Displaced dileptons 
only occur for our GMSB Higgsino models, in which case the long-lived particle is always 
neutral.) There is still some question about how this analysis would behave when additional 
tracks emerge from the same vertex, such as h ^ ZZ* ^ Z+/“-|-jets from Higgsino decay. 
However, this model is the only case that we study where such a question would arise, and 
the limits from this search are regardless not the most powerful for displaced Higgs decays. 


D. CMS Displaced Electron and Mnon 

The GMS displaced electron and muon search [T^] uses a highly minimalistic and inclusive 
analysis, simply demanding the presence of exactly one electron and one muon, each with 
signihcant 2D impact parameter up to 2 cm. The analysis is broken down into three exclusive 
bins in joint impact parameters, with varying degrees of background. Gonsequently, the 
statistical analysis is somewhat more involved, as discussed below. 

The electron and muon must each have pt > 25 GeV, \t]\ < 2.5, and transverse impact 
parameter between 0.02 cm and 2 cm. The last requirement allows the analysis to focus on 
the region of tracking parameters where the efficiencies are largely the same as for prompt 
tracks. The two leptons must also be well-isolated. In addition to a basic particle isolation 
requirement, which we form by tallying all hadrons (charged and neutral) within AR < 0.3 
and demanding a relative pt less than 10%, each lepton must also be isolated from jets with 
a Pt threshold of 15 GeV within AR < 0.5. The leptons themselves must also be separated 
by at least this distance and have opposite charges. 

The three exclusive signal regions consist of a high-displacement/lower-background region 
SR3 where both leptons’ impact parameters are above 0.1 cm, an intermediate region SR2 
where the event fails this criterion but still has impact parameters above 0.05 cm, and a low- 
displacement/higher-background SRI where the event fails both of these criteria but still 
has impact parameters above 0.02 cm. The observed (expected) event counts are, respec- 
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tively, 0 (0.051), 0 (1.01), and 19 (18.0). We have coded this three-bin statistical analysis 
as a toy Monte Carlo using a Poisson likelihood-ratio discriminator built from the central 
background predictions, and including the systematic uncertainties on the background (as¬ 
sumed Gaussian and uncorrelated) as perturbations on the simulated pseudo-experiments. 
This allows us to map out the 95% CLs boundaries in the three-dimensional space of signal 
bin counts. Our statistical analysis is not exactly the same as that performed by CMS, but 
it should furnish an adequate approximation. We have verihed this by reproducing sections 
of the leptonic RPV stop limits presented in the analysis note. Approximately speaking, 
high-lifetime signals that are concentrated in SR2 and/or SR3 have a limit A^(SR2)-|-A^(SR3) 
< 3, whereas low-lifetime signals that are concentrated in SRI have a limit A^(SRl) < 13.5. 

Other than the nontrivial statistical analysis, this particular displaced search was one 
of the simplest for us to implement, since it is insensitive to efficiency degradations and 
nontrivial geometries that occur for decays in the body of the detector. Also, due to the 
fact that the focus is on decays that occur before traversing the pixels, the issue of whether 
the displaced particle is charged does not appear. 


E. ATLAS Muon Spectrometer 

The ATLAS muon spectrometer search [20] (7 TeV, 2 fb“^) is focused on models where 
particle decay lengths are several meters, and have high probability of decaying outside 
of the HCAL. It uses a novel vertex-hnding algorithm [63] to identify the sprays of tracks 
from a displaced decay within the muon chambers. Events with two successfully identihed 
candidates are used in the analysis, with zero events expected and zero observed (again, 
setting an upper limit of about 3 signal events). 

To pass the analysis cuts, hrst of all both decays must occur within hducial regions of the 
muon spectrometer. The identihed decays must be well-isolated from tracks above 5 GeV 
out to Ai? < 0.4, and from jets above 15 GeV out to Ai? < 0.7. (Again, we do not consider 
displaced decays from charged particles, since the particle’s own track would veto it.) 

This analysis appears to hll an important niche at lifetimes intermediate between the 
tracker radius and beyond the outer detector radius for colored long-lived particles, and 
uniquely probes out to the largest possible decay distances for long-lived particles that lack 
charged states. Nonetheless, we will see that the search is ultimately not very powerful. 
There are several reasons for this: 

• Unlike all of the other analyses that we study here, it has not (yet) been performed at 
the full beam energy and luminosity. We will indicate how much this might help below 
by making a naive projection to 8 TeV, 20 fb“^, assuming identical signal efficiencies 
and zero background. 

• The ATLAS data acquisition becomes highly inefficient for particles traveling large 
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distances at sub-light speeds, and loses the signal when a time delay of ~ 7 ns has 
accrued relative to a light-speed particle. For example, a particle traveling at 0.7c and 
decaying at r = 5 m (at 77 = 0) accrues this much delay, and this velocity is already 
above the typical median for pair-produced heavy particles. The timing requirement 
therefore usually has an 0 ( 1 ) impact on our signal acceptance. 

• The analysis requires coincident behavior for both decays, and therefore pays every 
possible inefficiency factor twice. These include basic requirements such as neutral 
R-hadronization probabilities, the geometric constraints of the muon spectrometer, 
the vertexing efficiency, and the isolation requirements. The geometry in particular 
becomes a major factor when considering either very long or very short lifetimes, 
where it respectively leads to either a power suppression or an exponential suppression. 
Analyses that can rely on one candidate only need to be “lucky” once. 

Regarding the last point, it would be very interesting to see if this analysis could be 
run with a single-candidate option. On the one hand, this would result in much higher 
background rates (effectively 0(10 pb)) with the original reconstructions and cuts. On the 
other hand, many of the models that we consider here have quite appreciable cross sections, 
and because they have much higher masses than the baseline models that ATLAS studied, 
should lead to even more spectacular multi-track signatures in the muon spectrometer. 
Relaxing the isolation requirements somewhat, to allow more of the decay particles to point 
back to the HCAL, could also be benehcial. These HCAL signals would also be rather 
distinctive given that they would contribute mainly in the outermost layers. However, 
it is unclear whether ATLAS’s jet reconstruction requirements would anyway ignore such 
anomalous deposition patterns. 

It should be also noted that the information on the performance of the displaced vertex 
reconstruction from ATLAS’s papers is limited to a rather small set of new physics models, 
using the common benchmark scenario of a Higgs-like scalar decaying to a pair of displaced 
pseudoscalars. Only four mass points with fairly similar kinematics are studied, the most 
energetic decays coming from a 140 GeV scalar decaying into a 40 GeV pseudoscalar. It 
is therefore unclear how this very complex search would perform on, say, a 1 TeV RPV 
gluino decaying into three jets. Moreso than most of our other recasts, this one must be 
then viewed cautiously and somewhat conditionally. Still, because of the search’s limitations 
discussed above, in practice it does not end up probing masses beyond a few hundred GeV. 


F. ATLAS Low-EM Jets 

The ATLAS low-EM jets search [19] focuses on coincident jet-like signals conhned entirely 
to the HGAL, with stringent cuts on nearby EGAL and tracker activity. It is sensitive to pairs 
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of displaced decays within the HCAL volume. The analysis has a small but non-negligible 
background, leading to an upper limit of about 20 signal events. 

For a pair of jets to pass into the analysis, the leading (subleading) jet must have pr > 
60 GeV (40 GeV). Each jet should have no associated tracks with pt > I GeV within 
AR < 0.2, and should have at least approximately 16 (10^'^) times more energy recorded 
in the HGAL than in the EGAL. As discussed in Appendix our default model of these 
isolation requirements uses a combination of a flat efficiency factor and an over conservative 
veto on decays in the HGAL body that produce particles pointing back to the EGAL. An 
alternative, looser version removes the latter requirement, and the two extremes define an 
approximate error band for our modeling of this analysis. We also automatically veto events 
containing charged long-lived particles, which as usual would leave a track. ATLAS further 
makes an explicit cut of 5 ns on the signal timing delay relative to what would be deposited 
by a particle moving at light-speed. Given that the linear distances involved are about half 
as large as those for the muon spectrometer analysis, this ends up having a relatively less 
detrimental effect (though still potentially 0(1)) on the signal acceptance. Finally, there 
is a requirement that the event has < 50 GeV in order to reject non-collision events, 
though this can severely impact many of our models with invisible LSPs. 

Similar to the ATLAS muon spectrometer search, the requirement of coincident decays 
within the same detector system significantly limits the model reach. For our SUSY models, 
much of the candidate-by-candidate inefficiency comes from our somewhat ad hoc require¬ 
ment that no decay particles can point back to the EGAL, and we take results with and 
without this requirement to define our uncertainty band. Still, we conclude that this search 
is not very competitive, even to a luminosity-scaled muon spectrometer search. It seems as 
if it is simply too constrained by geometry. 

As in the previous subsection, it would potentially be useful if the search could be adapted 
for single-candidate acceptance instead of requiring decays with coincident properties. The 
fact that this would increase the background could be offset in several ways that should 
maintain high signal acceptance. First of all, most of the SUSY models that we consider 
would deposit far more energy than ATLAS’s benchmark models. A search that is broken 
down into different jet E^p bins could already improve sensitivity. Second, these very en¬ 
ergetic and wide-angle multibody decays might leave quite unusual 3D spatial and timing 
deposition patterns in the HGAL, possibly so unusual that even the relative EGAL deposi¬ 
tion requirement could be relaxed. In fact, we are again already giving ATLAS the benefit 
of the doubt by assuming that such unusual jets would pass the reconstructions of their 
analysis. But the limited mass reach of the search by itself limits the possible impact of this 
subtlety. 

It would also be very advantageous if the cut could be eliminated. For GMSB and 
mini-split models in particular, the cut is a major handicap. It can also contribute a subtle 
geometric problem in models without true ^t- Under the zeroth-order assumption that a 
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particle decaying within the HCAL has all of its energy absorbed at one point, the particle’s 
reconstructed momentum vector is effectively rescaled by its lab-frame energy. For sparticle 
pair production that is not exactly back-to-back in 3D space, true transverse momentum 
balances, but the energy-scaled transverse momentum need not. Again, this is mainly an 
issue for models with large masses, beyond the sensitivity of the original search. But it could 
become problematic if the search were to be extended. 


G. ATLAS Displaced Muon Plus Tracks 

The ATLAS muon plus tracks search [T3] uses the inner tracker to reconstruct highly 
displaced vertices containing at least one muon.® The basic search, which counts the number 
of events containing at least one such vertex, is background-free. However, a looser version of 
the search relaxes the demand that the muon is matched to the displaced vertex candidate. 
This version is also background-free, and has improved signal acceptance, especially for 
models with a mixture of low-multiplicity leptonic decays and high-multiplicity hadronic 
decays (i.e., our GMSB Higgsino and stop models). We utilize this version for our recasts. 

Vertices are reconstructed from displaced tracks originating within the inner region of the 
tracker, r < 18 cm and \z\ < 30 cm. The tracks used toward the vertex reconstruction must 
have pt > 1 GeV and transverse impact parameter greater than 2 mm. A good candidate 
vertex must be reconstructed within the hducial volume, have at least hve associated tracks, 
a track-mass greater than 10 GeV. The muon selection is pt > 55 GeV and |p| < 1.07, and 
transverse impact parameter greater than 1.5 mm. 

Unlike the other ATLAS searches, this one requires neither very tight activity cuts nor 
coincident behavior for the decays. Gonsequently, it is much more powerful within its realm 
of applicability. Again, we conservatively ignore displaced particles produced in the decays 
of charged R-hadrons. 

Because this is a tracker-based analysis, similar to the GMS displaced dijets, the possi¬ 
bility also exists for subtleties when heavy flavor secondary decays occur after the displaced 
decay. However, ATLAS states explicitly that vertices less than 1 mm from one another are 
merged. While the exact behavior for heavy flavor hnal-states is not given by ATLAS, we 
assume that this merging procedure effectively makes them insensitive to this issue. 


III. MODELS AND LIMITS 

We now describe three of the well-known scenarios that lead to displaced sparticle decays, 
and present our new limits on several simplified models within those scenarios. In all three 

® We do not consider the superceded versions of this search at lower luminosity and beam energy [BH [BS] . 
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scenarios, we consider models with LSP or NLSP glninos, which are common search targets 
dne to their enormous pair production cross sections. For gauge mediation, we also consider 
the closely-related case of light-squark NLSPs. Otherwise, we focus on either stop pair 
production or Higgsino pair production, as both particles are expected to be below 1 TeV 
in a truly natural SUSY theory In most of what follows in this section, we restrict 

ourselves to presenting the basic search results, and reserve commentary on implications for 
Section HVl 

We generate event samples for most models with PythiaS [58], making extensive use 
of that program’s R-hadronization capabilities. Final-state particles from each lowest-level 
R-hadron or Higgsino decay are subsequently displaced before detector simulation and event 
reconstructions. For some of the colored production models where multibody decay kine¬ 
matics can be important, we have generated events in MadGraphS |HS]- In this case, the 
PythiaS R-hadronization routine does not work because the R-odd colored particle has al¬ 
ready decayed, and its daughters given color connections to other parts of the event. There, 
we identify final-state hadrons as descendants of the long-lived colored particles in an ap¬ 
proximate way: each hadron is associated to the closest quark/gluon in AR as viewed the 
end of the parton shower, and the ancestry of that quark/gluon is traced by proxy.® 

All pair production cross sections are normalized to their NLO-I-NLL predictions, includ¬ 
ing colored production through pure QCD [67] and electroweak Higgsino production [68] . 
These predictions are all conveniently tabulated for 7 and 8 TeV by the LHC SUSY Cross 
Section Working Group [Hn]- The Higgsino predictions assume nearly mass-degenerate Hig¬ 
gsino states with small mixings into the gauginos (assuming Mi 2 ~ 1 TeV). Those cross 
sections are only provided up to 410 GeV, but we assume a flat K-factor for higher masses. 

The generated particle-level events are passed through the detector simulations described 
in the Appendix and then subjected to the various analysis cuts described in Section |TT| 
The main output is a set of experimental acceptances for each individual simplified model 
(either per-decay or per-event), scanned over mass and lifetime. As an illustrative example, 
we provide our acceptances for the RPV i ^ ds model described below, passed through 
the GMS displaced dijet analysis. As can be seen, this specific analysis is most efficient for 
masses greater than a few hundred GeV and lifetimes at the 0(10 cm) scale. 

When constructing limits, for most of our recast experimental searches we provide a 
very rough error band on our predicted exclusion regions by varying these signal acceptance 
estimates up and down by a factor of 1.5 (not allowing those rescaled acceptances to exceed 

^ We save explicit investigation of a left-handed sbottom (N)LSP for future work. 

® Using the simple 2-body t —)■ didj as a cross-check, we find that the invariant mass of each stop is 
reconstructed to within about 10% and without bias, and that most search acceptances are only mildly 
sensitive to to the ordering of decay and hadronization. The largest effects are on the tracker-based 
searches, with the efficiencies dropping by 10-15% for the sample decayed before hadronization, and 
therefore yielding conservative limits. 
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FIG. 1: Our central predictions for the CMS per-decay displaced dijet acceptances for t ds, using 


The different colored lines are steps of 100 GeV in mass. 


the detector model described in Appendix and following the analysis described in Section IIB 


unity). This gives some indication of sensitivity to possible recasting errors. There are only 
two specihc searches where we do not follow this protocol. The hrst such search is for stable 
charged particles, for which we do not explicitly include an error band. Our modeling here is 
fairly basic and conservative, and the acceptance anyway turns off exponentially fast at low 
lifetimes. We have also recast CMS’s over-conservative “charge stripped” limits, estimated 
in a scenario where interactions in the calorimeters always strip off the R-hadron charges. 
The second search is for low-EM jets at ATLAS, where we have opted to instead dehne 
over-conservative and under-conservative treatments of the isolation against EM calorimeter 
activity, which we cannot reliably model. Here, we require either that no decay particles 
point back to the ECAL, or do not place any explicit isolation criterion (though in both 
cases we employ a flat 0(1) reconstruction efficiency factor given in Appendix F). 
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A. Baryonic R-Parity Violation 


One of the simplest extensions to the MSSM is the introduction of R-parity-violating 
Yukawa superpotential couplings and/or a /i-term between the lepton doublet and down- 
type Higgs doublet superhelds [3]. R-parity violation may also be introduced in the soft 
SUSY-breaking potential, or in the Kahler potential (HD]- These are all typically set with 
zero coefficients in order to enforce R-parity. R-parity trivially prevents dimension-four 
proton decay and stabilizes the LSP, providing a possible dark matter candidate. However, 
these virtues are hardly strict requirements of the MSSM. Proton decay can be prevented 
by alternative stabilizing symmetries, which in any case may be required given the existence 
of potentially dangerous R-even operators at dimension-five mi- Dark matter could easily 
arise from a different particle sector instead of the MSSM neutralino. 

Violation of R-parity can lead to radical changes in collider phenomenology, depending 
sensitively on which operators are activated, on the magnitude and flavor structure of those 
operators, and on the identity of the LSP, which is no longer stable and no longer needs 
to be electromagnetically neutral. Broad ranges of coupling strengths allow for the LSP 
to decay at displaced locations within the LHC detector volumes. For example, for two- 
body sfermion decays into light SM fermions, mediated by one of the R-parity-violating 
Yukawas, a dimensionless coupling of would yield a substantial population 

of measurably-displaced decays. 

Proton stabilization allows the active RPV operators to violate either lepton number 
or baryon number, but not both. This partitions the RPV scenarios into two mutually- 
exclusive classes, which we can call leptonic RPV and baryonic RPV. The only explicit 
RPV displaced searches so far at the LHC have assumed leptonic RPV [m [151 113 HH], 
capitalizing on the presence of leptons in the final state to help with triggering and with 
controlling backgrounds. Here we pursue the completely untested case of baryonic RPV. 
(However, we will put the leptonic RPV searches to an alternative use in the next subsection, 
applying them to gauge mediation models.) 

We begin with the case of stop LSP. Minimalistically, the decay would be mediated 
by the usual Yukawa superpotential interaction of right-handed chiral quark/squark fields, 
X'-ji^UiDjDk {ijk are flavor indices). This is the only baryon-number-violating interaction 
in the MSSM that respects SUSY at dimension-four. The required antisymmetrization over 
color indices requires a commensurate antisymmetrization over down-type flavor indices, 
leading to allowed decays i —>■ ds, db, and sb. Recently, it has also been observed that 
stop decays may proceed through a different combination of chiral quark/squark fields, via 
dimension-five operator QiQjD\. in the Kahler potential [HO]- The resulting component-field 
operators allow for a decay i —)• bb, and indeed this is generally preferred since the decay 
amplitudes are chirally-suppressed (analogous to pion decay). Prompt decays i —)■ djdk 
have only just begun to be probed by the LHC, in the mass region 200-400 GeV m- It 
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FIG. 2: Recast constraints on displaced t —>• djdk via baryonic RPV. Colored bands indicate 
acceptance variations up/down by 1.5. The dot-dashed lines indicate contours of ^^ 12 , assumed to 
be the only contributing RPV coupling. Prompt limits (dark gray) are from [70], and low-mass 
search projections (light gray) are from [^. They are conservatively cut off at 1 mm. 


has been estimated that a search based on jet substrncture could also push down into the 
lower-mass region currently not covered [7T]. (For longer-term projections, also see [7T] as 
well as [72].) The only other available limits are when the stops are detector-stable, the 
strongest (~ 900 GeV) coming from the CMS and ATLAS charged R-hadron searches [318]. 

Figs. 1^ and show the regions of mass and lifetime for i —>■ djdk that have now been 
excluded according to our recasts, taking the two extreme cases of only light-flavor decays 
and only bb decays. The sensitivity is dominated by the charged R-hadron and displaced 
dijet searches, a pattern that will recur often in our colored sparticle limits. For both 
models there is nearly complete coverage out to almost 1 TeV, with a notable weak-spot at 
cr ~ 10 m and of course much weaker limits for displacements < mm. This weakening at 
low lifetimes is more pronounced for the bb decays, partially because the CMS dijet search 
is intrinsically less efficient for heavy flavor decays due to the somewhat lower particle track 
multiplicities, but also because of the conservative choice in our modeling of displaced vertex 
reconstruction for 6-jets, discussed in Section IIB At lower lifetimes, we have also indicated 
the existing and projected prompt limits, applying a conservative sensitivity cutoff at 1 mm. 
(There should still be sensitivity from prompt searches for longer lifetimes, but we do not 
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FIG. 3: Recast constraints on displaced t —> 66 in the Dynamical RPV framework. Colored 
bands indicate acceptance variations up/down by 1.5. The dot-dashed lines indicate contours of 
r/ggg/M, assumcd to be the only contributing RPV coupling. It arises from the Kahler potential 
operator M)Q'j,Q‘iD\ + (h.c.). Prompt limits (dark gray) are from [70] (neglecting a possible 
improvement in the limits due to the higher b-jet multiplicity in the DRPV model), and low-mass 
search projections (light gray) are from [71]. They are conservatively cut off at 1 mm. 


have enough information to reliably model this.) Combining these three searches, unbroken 
coverage is achieved for all lifetimes for masses where the prompt searches are sensitive. 
Indeed, for stop masses up to a few hundred GeV, the CMS dijet search alone spans 6-7 
orders of magnitude in lifetime. This amazing performance capitalizes heavily on the fact 
that millions of stop pairs would have been produced at such small masses, with sizable 
enough kinematic tails to pass the jet Ht and pt cuts, and enough remaining rate to catch 
anomalously early or late decays from models at the edges of the exclusion region. Two 
other displaced searches, ATLAS muon chamber (including our naive 8 TeV projection) and 


ATLAS low-EM jets, are much less competitive for the reasons discussed in Sections IIE 


and jll F , though they do offer useful complementarity in that their limits are derived from 
completely different detector systems. Finally, we point out a sizable region in the bb decay 
case that is also covered by the ATLAS /i-|-tracks search, from a small population of events 
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FIG. 4: Recast constraints on displaced g —)• jjj via baryonic RPV. Colored bands indicate 
acceptance variations up/down by 1.5. The dot-dashed lines indicate contours of We 

have parametrized the decay assuming that one species of off-shell RH squark dominates, and splits 
into quarks via a single A" ;, coupling. All final-state quarks are also assumed to be from the first 
two generations. Prompt limits (gray) are derived from [73]. They are conservatively cut off at 
1 mm. 

where one of the bottom decays produces a hard muon.® 

The next model that we consider is gluino LSP. Considering only traditional superpoten¬ 
tial RPV, the gluino decays by hrst transitioning into a virtual squark and a corresponding 
real quark. The virtual squark then splits to two quarks through the UDD operator. The 
full 3-body decay is ^ jjj. There are again many options for flavor structure, which 
may be engineered both at the level of the A" ;, couplings and the squark mass spectrum. 
Here, we simply assume decays into light flavors, though decays involving 6-quarks could 
again be subjected to weaker limits at low lifetimes, and decays involving t-quarks would 
also receive constraints from the displaced searches involving leptons. Otherwise, we expect 
fairly similar coverage. Of course, branching ratios into top also suffer additional phase 
space suppression. 

Fig. a shows our estimated exclusions for g —)■ jjj. The qualitative features are quite 

® The muon in this search is triggered from the standalone muon spectrometer, and is not explicitly required 
to be isolated. 
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FIG. 5: Recast constraints on displaced —)• jjj via baryonic RPV. Colored bands indicate 

acceptance variations up/down by 1.5. The dot-dashed lines indicate contours of rriq /We 
have parametrized the decay assuming that one species of off-shell RH squark dominates, coupling 
to the Higgsino according to its up-type or down-type Yukawa Yg oc nig/{v sin 13) or mg/(u cos/3), 
and splits into quarks via a single coupling. All final-state quarks are also assumed to be from 
the hrst two generations. 

similar to the RPV stop decays, though the much higher cross sections yield a signihcantly 
extended mass reach for all searches. CMS dijets in particular reaches above 1.5 TeV, close 
to the production limit of 0(1) event in the entire run, and exceeding the mass reach of the 
stable R-hadron search by several hundred GeV. Notably, the displaced trijet conhguration 
is very efficiently picked up by the CMS dijets search, which was designed for a very different 
signal. The weak spot at 10 m is still apparent, but much less pronounced since the CMS 
dijet search nearly matches the HSCP search sensitivity at that lifetime. It is also interesting 
to supplement with the limits from prompt searches 1131 EH, which are similar for purely 
light-flavor decays and decays containing 6-quarks. Again applying an ad hoc 1 mm cutoff 
on the lifetime sensitivity of the prompt searches, there is currently unbroken coverage for 
all possible lifetimes for masses potentially as high as 900 GeV. 

The last baryonic RPV example model that we consider is a Higgsino multiplet “co- 
LSP.” The four Higgsino states are assumed to be only mildly mixed into heavier electroweak 
gauginos, and the multiplet split by 0(10 GeV) or less. The heavier Higgsinos undergo a soft 
but prompt cascade via virtual gauge boson emission into the lightest, neutral Higgsino. The 
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displaced decay of this lightest Higgsino proceeds in almost exact analogy with the gluino 
decay, though now the virtual squark is accessed via a super-Yukawa coupling instead of a 
super-QCD coupling. Again, the flavor structure of the decay can be nontrivial, but as a 
first pass we simply assume that the final quarks are all light. Decays involving tops deserve 
a dedicated investigation, especially in the context of a natural theory, though we anticipate 
fairly similar limits. 

Fig. I shows our estimated exclusions for H —)■ jjj. The qualitative picture from the 
displaced decay searches remains similar, though with reduced mass reach due to the smaller 
production cross sections. Unlike the preceding stop and gluino examples, there are no 
explicit limits on —)■ jjj in either prompt or stable charged particle searches. For the 
prompt case, we can compare to ^ jjj searches Hal [71]. Using simple cross section scaling 
suggests that that promptly decaying RPV Higgsinos are genuinely unconstrained, since the 
Higgsino cross section is roughly 500 times smaller at a given mass. (A more aggressive 
dedicated prompt search could be useful, though would be highly challenging.) For the 
stable case, the LSP here is generically neutral, and hence does not leave a track. Therefore, 
our reported direct Higgsino production limits here are the first for any lifetime. 

While we have only studied a small sample of possible spectra, these results clearly 
illustrate the power of the LHC in probing baryonic RPV in general via displaced decays. 
An obvious extension of our observations would be an application to a broader class of 
flavor assumptions, though as indicated we do not expect radically different sensitivity. 
The remaining extensions would be to consider different LSPs, and perhaps more model- 
dependent scenarios where the LSP is created in cascades in addition to prompt production. 
An LSP squark could represent a rather trivial example, since the production and decay could 
be very similar to the LSP stop. However, decays into the top quark could also open up, 
and the effective production cross section could also be highly enhanced by the multiplicity 
of nearby squark states (cascading promptly into one another) and/or by gluino exchanges. 
Direct LSP slepton production represents a qualitatively different direction, wherein a 4- 
body decay to Ijjj or ujjj (via virtual electroweakino and squark) might dominate, even 
for much larger values of the A" ;.. Finally, production of different electroweakinos, such as 
a mostly-bino or mostly-wino, could be considered. In fact, the latter has recently been 
investigated in |3l|, and similarly hnds very high mass reach using the CMS displaced dijets 
search. With the generality and power of the HSCP and dijets searches, the main missing 
pieces in covering the mass-lifetime plane for these varied models would be prompt and b- 
tagged searches (possibly recast from other models) and more aggressive muon-chamber and 
calorimeter searches, especially for the cases without long-lived charged states. Additional 
studies within the framework of dynamical RPV could also be interesting, since this allows 
for additional flavor and chirality structures in the decays. 
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B. General Gauge Mediation 


Traditional forms of gauge mediation assume fairly minimal messenger sector construc¬ 
tions, and consequently relatively predictive patterns of sparticle masses |1]. For example, 
relationships like Mi.M 2 '.M^ = ai:a 2 '-(y 3 for the three gauginos favor a bino-like NLSP and a 
much heavier gluino. A much more general perspective has been advocated in izg, acknowl¬ 
edging the full range of possible MSSM spectra derivable from arbitrary messenger sectors, 
perturbative or not. Practically speaking, the individual sfermion and gaugino soft masses 
become almost freely-adjustable, up to two sfermion sum rules and flavor universality, as well 
as vanishing A-terms at the mediation scale. Even more general model frameworks allow 
for the possibility of flavor nonuniversality effects, either by mixing into a supersymmetric 
composite sector [76HT9], by gauging flavor symmetries [HQ], or by introducing large A-terms 
through non-minimal interactions between MSSM and messenger fields [8T] . 

Such freedom of model-building in gauge mediation allows for a number of intriguing 
options for LHC phenomenology. Practically any superparticle can be made the NLSP, which 
then decays into its SM-partner and a light gravitino (i.e., Goldstino) at a rate controlled 
only by the fundamental SUSY-breaking mass scale, a/F. Roughly speaking, when the SM 
partner is light. 


CT ~ 0.3 mm 


100 GeV 


m 


5 


Vf \ 

100 TeV j ’ 


( 1 ) 


implying displaced decays at the 0(mm - 10 m) scale for y/F in the range of a few-hundred 
to a few-thousand TeV. With the traditional option of the bino-like neutralino as the NLSP, 
the decays are dominated by photon and gravitino. Gonsequently, the experimental effort in 
displaced GMSB searches has focused on signals of displaced/delayed photons and missing 
energy In ^m\, it was pointed out that an NLSP neutralino with larger wino or 

Higgsino fraction would also yield displaced W, Z, and Higgs. The possibility of displaced 
NLSP stops was emphasized in jHU |83l El], especially for stop masses near or below the top 
quark mass. The simplest remaining displaced NLSP options would be slepton decaying to 
lepton, squarks decaying to (non-top) quarks, and gluino decaying to gluon. None of these 
other possibilities have been searched for in the case of displaced decays, though a number of 
searches have been performed for prompt decays and for the collider-stable cases. However, 


These are also necessarily the most model-dependent of the available signals, since a simplified model 
containing only the bino would have vanishing tree-level pair production cross section. The most powerful 
existing search m mainly capitalizes on pair production of mostly-winos cascading down into a mostly- 
bino, relying on the relationship M 2 /Mi ~ a 2 /oii. 

An NLSP sneutrino from light doublets would decay fully invisibly, making the displacement irrelevant 
for its experimental signatures. The distinctive phenomenology of such a scenario has been studied in |85P 

E7]. 
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most NLSP possibilities are actually already under tight constraint, as we will see, again 
sometimes with coverage over the full range of possible lifetimes. 

A notable exception is any variation on the slepton NLSP, such as the standard mostly- 
tr when the sleptons are degenerate up to Yukawa effects, or possibly eu/jlR if the staus 
can be made heavier |2S]- These would be largely unconstrained by the displaced vertex 
searches due to the low track multiplicities and vertex masses, and searches for displaced 
activity in the calorimeters or muon chambers would usually fail to pick up the signal, for 
example because the associated slepton track would cause isolation failure. CMS’s displaced 
e+fi search na should yield some sensitivity for leptonic tan decays. Though we reserve for 
future work a more comprehensive study of the status of displaced slepton NLSPs in GMSB, 
we anticipate that planned searches for “kinked track” topologies will be important to more 
fully cover the parameter space, and that existing disappearing track searches might also 
provide some sensitivity. 

The remaining NLSP options that yield a single SM final-state particle, without passing 
through an intermediate heavy SM decay (top or electroweak), are the non-top squarks and 
the gluino. As in the baryonic RPV models, presently the only potentially applicable tracker- 
based search is the CMS displaced dijets, but the nominal number of partons in the decay is 
not two. Here, in order for the displaced dijet search to be sensitive, the decay must undergo 
a hard enough final-state bremsstrahlung to create a second jet. Given the large strong 
production cross sections, this is an affordable penalty: of order (a^/vr) log(m/(60 GeV)), 
times color factors, with the sparticle mass and jet pt cut appearing inside the logarithm. 
Quantitatively, the chance to radiate a second jet is roughly at the 1-10% level. Of course, a 
search explicitly geared toward the one-jet topology could be more efficient, and the displaced 
dijet trigger would already capture this signal in cases where both decays occur within the 
inner 60 cm radius of the tracker and with at least 0(mm) displacement. However, at the 
very high and very low lifetime ranges, the inefficiency induced by requiring an extra jet 
may be less than the inefficiency that would be induced by forcing both decays to occur at 
improbably short or long proper times. A more fruitful option for future analyses could be 
to exploit the traditional jet. Hr, and triggers, and apply an offline search for individual 
displaced jets.^^ 

Proper simulation of the decays for the existing displaced dijet search requires some level 
of matrix element matching. This is performed automatically by the PythiaS shower in the 
case of g —>■ with a massless neutralino LSP, while the desired decay q —)■ qG is not 

matched. However, we have observed essentially identical rates and kinematics for extra jet 
production in q decays between explicit MadGraphS 2- and 3-body decay simulations with 
neutralino and gravitino LSPs, and close agreement with PythiaS’s matched predictions for 
the first shower emission. We therefore feel confident using the massless neutralino LSP as a 


We thank Joshua Hardenbrook for emphasizing this possibility. 
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q-*qG (GMSB) 



ifiq (GeV) 


FIG. 6: Recast constraints on displaced q —)• qG in general GMSB, conservatively assuming con¬ 
tributions from only (Ir and sr. Colored bands indicate acceptance variations up/down by 1.5. 
The dot-dashed lines indicate contours of the SUSY-breaking scale ^/F. Prompt limits (gray) are 
derived from [89]. They are conservatively cut off at 1 mm. 

proxy for the gravitino LSP for squark decays in PythiaS. For the gluino, such an analogous 
decay to neutralino does not exist at tree level, is not part of the MadGraphS MSSM model, 
and would not obviously be matched if forced to proceed in PythiaS. Instead, we compare 
the unmatched PythiaS predictions for its hrst shower emission to MadGraphS, both with 
gravitino LSP. We again hnd similar decay kinematic distributions, with PythiaS predict¬ 
ing a somewhat slower falloff out to AR{j,j) ~ n. But the major difference is in the total 
emission rate, which PythiaS over-estimates by a factor of about 1.8. To approximately 
compensate for this, we rescale the individual vertex reconstruction efficiencies by 1/2. It 
should be understood that 0(10%) modeling uncertainties on the displaced dijet reconstruc¬ 
tion efficiencies for GMSB gluinos should likely be applied, though we anyway effectively 
absorb this into our ad hoc systematic variations. 

Starting with the squark NLSP, we display the results in Fig. We conservatively 
assume just two degenerate species, cIr and sr. This is a technical possibility if the SU{3) 
contributions to the sfermion masses are small, the SU{2) contributions are large, and the 
third-generation squarks receive additional mass contributions. The exclusions are similar to 
those of the RPV stops (Figs. [^and[^, although now with much stronger prompt jets-|-.^T 
searches. Unbroken coverage over lifetime is achieved up to about 450-550 GeV, limited by 
the crossover between the HSGP and displaced dijet searches. 
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g -* g G (GMSB) 
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FIG. 7: Recast constraints on displaced g — )• gG in general GMSB. Golored bands indicate accep¬ 
tance variations up/down by 1.5. The dot-dashed lines indicate contours of the SUSY-breaking 
scale y/F. Prompt limits (gray) are derived from the q —)• qG of [89j. They are conservatively cut 
off at 1 mm. 

Next we consider the gluino NLSP in Fig. Comparing to the RPV results for g —)■ jjj 
in Fig. 1^ we observe much weaker displaced decay limits and much stronger prompt decay 
limits. The former is due to the requirement of additional hard radiation in the decay to 
pass the CMS displaced dijet reconstruction. The latter is due to the much more distinctive 
jets-|-.^T signature. Most of the model space below 1200 GeV is covered, with expected weak 
spots at 0(mm) and 0(1-10 m), though much of the surviving space at smaller lifetimes 
would likely be probed by a more detailed iets+^x recast, as in [10]. Full coverage over all 
lifetimes is only achieved for masses below 800 GeV. 

We now move on to the naturalness-motivated options, starting with the NLSP stop 
in Fig. 1^ We consider stops of any mass above 100 GeV, including a range of masses 
below rrit and through the compressed region where In these regions, the 

decays are dominantly 3-body t —)• WhG, with a large fraction of energy going into the 
effectively derivatively-coupled gravitino/Goldstino. Also, in addition to the by-now familiar 
searches that have appeared in all of our recasts above, the semileptonic decays of the stop 
open up sensitivity in the ATLAS /i-|-tracks search and the GMS displaced e-f/i search. 
While the sensitivity regions for these searches are fully contained by GMS displaced dijets. 


Because of the smallness of the ttG coupling, top quark decay into a light stop and G would be rare. 
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t-»tWG (GMSB) 



FIG. 8: Recast constraints on displaced i —> t^*^G in general GMSB. Golored bands indicate ac¬ 
ceptance variations up/down by 1.5. The dot-dashed lines indicate contours of the SUSY-breaking 
scale \/T. Prompt limits (gray) are derived from |9n[l91j . They are conservatively cut off at 1 mm. 

corroborating coverage is provided by the leptonic searches over much of the excluded region. 
Adding in the prompt searches [901191], which likely give unbroken coverage between 100 GeV 
and 670 GeV,^"^ we infer that GMSB stops of any lifetime are excluded below about 500 GeV. 
For lifetimes at the cm-scale, exclusions extend beyond 700 GeV, and, as noted before, out 
to about 900 GeV for lifetimes longer than 0(10 m). (For other estimates of displaced stop 
exclusions in GMSB, see [35].) 

Finally for GMSB, we consider Higgsino multiplet co-NLSPs. As in the RPV case 
above, we assume that all Higgsino states are nearby to one another (split by no more 
than 0(10 GeV)), with heavier states decaying promptly. The lightest Higgsino will pref- 

The prompt searches face some subtleties. On the one hand, for stop masses well above rrit, existing 
searches for t with massless neutralino should offer identical coverage. On the other hand, the decay 

kinematics near or below m* can be significantly different than the corresponding decays to neutralinos. 
The expectation is that the GMSB limits there should be much stronger than the nominal limits, and 
not subject to the usual sensitivity gap with a compressed spectrum [S^j. The major exception is a 
mostly-^L stop, for which spin effects would reshape the .^T-sensitive distributions and weaken the limits 
in searches with semileptonic decays. Without recasting the most recent searches, it is not possible to 
precisely delineate this gap, though the results of ^ suggest that it may be several lO’s of GeV wide, 
centered in the vicinity of 200 GeV. Since the displaced searches are not designed to cut on .^T-sensitive 
tails, we do not expect such spin effects to be significant there. 
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a) H-*ZG (GMSB) 
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c) H->Z/hG, large tanp (GMSB) 
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FIG. 9: Recast constraints on displaced decays via GMSB: a) pure —)■ ZG, b) pure —)• 

hG, c) large-tanGolored bands indicate acceptance variations up/down by 1.5. The dot-dashed 
lines indicate contours of the SUSY-breaking scale ^/F. Prompt limits (gray) are derived from [93]. 
They are conservatively cut off at 1 mm. 
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erentially decay to some mixture oi ZG and hG, with 7 G suppressed. The lifetime and 
branching fractions of the lightest Higgsino exhibit simple behavior if mixings with the bino 
and wino are small, and the scalar Higgs sector is close to the “decoupling limit.” For 
instance, when tan/? = 1 , the lightest Higgsino coupling to either ZG or hG vanishes, de¬ 
pending on the relative signs of the /i parameter and Mi_ 2 . For tan/? ^ 1, the hG and ZG 
decay modes have similar partial widths if ^^h- 

We present the limits in these three extreme cases in Fig. a) pure H —)• ZG, b) 
pure H —)■ hG, and c) large-tan /?. Though the HSCP searches again do not apply, the 
GMSB Higgsino brings into play all of our other displaced decay recasts, now including as 
well the CMS displaced dilepton search. For decays that include direct Z bosons, this last 
search can be seen to play a major role, competing significantly with and even beating the 
CMS dijet search. This is especially obvious at very low and very high lifetimes, where 
the dilepton search benefits from both lower track impact parameter thresholds and looser 
Pt reconstruction thresholds. Similarly, leptonic decays play a major role in the prompt 
searches [HSj, with maximal sensitivity for pure H —)■ ZG, capitalizing in part on 3- and 
4-lepton channels. Taking the large-tan (3 case as a baseline example, rnp = 100 CeV is now 
covered from zero lifetime up to cr ~ 100 m, and mp = 300 CeV is covered up to cr ~ 
few m. The highest mass reach is for cr ~ 10 cm, extending out to about 600 CeV. For the 
pure H —)■ hG case, our displaced search recasts represent the only extant limits, as was the 
case for the RPV Higgsinos. 

We have touched upon most of the possible displaced NLSP decays in gauge mediation. 
The last obvious remaining option, which we now briefly discuss for completeness, would 
be wino co-NLSPs. In some ways, this should overlap significantly with our H results, 
but there are some notable differences. First, wino cross sections are about two times 
larger. Second, the has a significant branching fraction into photons, hence subjecting 
it to the displaced/delayed photon searches. Third, when the wino multiplet is somewhat 
separated in mass from the bino and Higgsino, the mass splitting between charged and 
neutral states is protected at leading order in the mass mixing by an accidental custodial 
symmetry, with the hrst nontrivial mixing contributions often comparable to or smaller 
than the electroweak radiative mass splitting of ~ 170 MeV. This famously leads to the 
long-lived decays hF^ —)■ , with cr ^ 5 cm, searched for in |T3l [55]. In such a case, 

there can be nontrivial competition between the above decay and —)■ W^G. There can 

also be peculiar cases with an initial stage decay —>■ , leaving a disappearing track, 

followed by a secondary displaced decay IV° —)■ {pf/Z)G. Finally, there are also some corners 
of parameter space with m{W^) < m{W^) due to chargino and neutralino mixings [94] . 
causing every event to contain two displaced IVs. The signatures would be much more 
similar to those of 77 —)• ZG, though missing the displaced dileptons and, if cr > 0(cm), 
containing a track or track stub pointing to the displaced decay (similar to the slepton 
NLSPs). All together, the potentially rich displaced phenomenology of wino co-NLSPs in 
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gauge mediation clearly deserves a more detailed investigation, and would bring together a 
surprisingly varied set of displaced search results. 


C. Mini-Split SUSY 

The last model framework that we consider is mini-split supersymmetry, where the scalars 
of the MSSM (excepting the SM-like Higgs boson) are all raised to the 1000 TeV scale [HI 
ii 1^ . This scale could represent a sweet-spot between masses that are high enough to 
avoid flavor constraints with arbitrary sfermion mass matrix structure, but low enough to 
provide a 125 GeV Higgs from the stop loop corrections. The separation between MSSM 
scalars and fermions can arise automatically in several SUSY-breaking mediation scenarios 
(surveyed in [5] and discussed on general terms in |6]). However, while some of the important 
virtues of SUSY such as gauge coupling unihcation and dark matter can be preserved, 
the original motivation of naturalness is partially abandoned. The apparently hnely-tuned 
Higgs mass might nonetheless be viewed as a byproduct of anthropic selection bias in the 
multiverse, in some ways similar to the unnaturally small cosmological constant [95], or as 
a compromise against much larger “tunings” within the available “model space” of broken 
SUSY theories |6|. 

While gauginos need not be present at any particular mass scale, the WIMP miracle 
is suggestive of TeV-scale masses, potentially within reach of the LHC. One of the most 
interesting targets is the gluino, since the flow of R-parity in its decay must pass through 
the heavy squarks, leading to suppressed matrix elements and extended lifetimes. Two types 
of LHC searches so far have directly targeted this signal: searches for R-hadrons stopped in 
the calorimeters and decaying out-of-time with respect to collisions [HI |2l| , and an ATLAS 
reinterpretation of its prompt jets+.^T searches using models with displaced decays [Tn] . 
Neither of these are optimally sensitive, though the former strategy has the added beneht of 
permitting a lifetime measurement if a positive signal is observed, and the latter strategy can 
be carried out quickly with no changes to the event reconstruction and selection software. 
In addition, searches for anomalous tracks from collider-stable R-hadrons, which we have 
discussed above in the contexts of both RPV and GMSB, continue to apply. 

Here we put all of these long-lived gluino searches into context for a couple of specihc 
assumptions for the decay kinematics. In almost full generality, the gluino may decay into 
any flavor-combination of quark-antiquark pairs plus a R, W, or ff. The exact admixture 
of these decays is highly model-dependent. Since the decay rate through any given off- 
shell squark channel scales like 1/m^, the lightest squark eigenstate would dominate if there 
is a somewhat spread-out scalar mass spectrum. Yukawa effects in the renormalization 
group may favor light third generation squarks, suggesting dominant decays g —>■ tix^ or 
g — )■ thx~+C.C. [6], though flavor-anarchic soft masses could instead lead to lighter hrst- or 
second-generation squarks. Loop-induced decays g —)■ gx^ might also become important [96] . 
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a) g -> q q B, m(B) = 0 (mini-split) 



Mg (GeV) 


b) g -> q q B, m(B) = m(g) - 100 GeV (mini-split) 



nig (GeV) 

FIG. 10: Recast constraints on displaced g —)■ qqB in mini-split SUSY: a) = 0, h) = 
nig — 100 GeV. Colored bands indicate acceptance variations up/down by 1.5. The dot-dashed 
lines indicate the intermediate squark mass, assuming that either dji or sr dominates the decay. 
Prompt limits (gray) are derived from [89(197]. They are conservatively cut off at 1 mm. Additional 
displaced limits come from stopped R-hadron searches [2Tj and ATLAS’s recast prompt limits [10] . 
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though again depending in detail on the squark mass spectrum, as well as on the gaugino 
mass spectrum. For our initial study here, we simply assume 100% 3-body branching fraction 
g —)■ qqB for q = udcs. The only free parameters are then the gluino and neutralino masses, 
as well as the gluino lifetime set by the squark mass scale [6], 


cr ~ 10 ^ m 


rrig /TeV\^ 
PeV/ \ rrig ) 


( 2 ) 


We reserve a more general survey of displaced mini-split phenomenology for future work. 

Fig. 10 shows our results, choosing either m^ = 0orm^ = rrig — 100 GeV. We hnd 
once again that, amongst explicit displaced decay searches, CMS displaced dijets offers 
superior sensitivity. For the light neutralino case, exclusions extend out as far as 1400 GeV, 
and for the heavy neutralino out to 800 GeV. The dedicated mini-split gluino searches, 
based on stopped R-hadron decays dSl IM] and recasts of prompt searches HD], do not 
tend to be competitive with this search combined with the HSCP searches. Interestingly, 
the ATLAS muon chamber and low-EM jets searches, which have tended to give universally 
weaker coverage, potentially offer more stable sensitivity as the visible activity gets squeezed 
out when ~ rrig. This owes largely to their focus on lower-mass signatures, which 
for compressed SUSY spectra becomes a major virtue. However, the overlap with HSCP 
coverage remains substantial. 


IV. DISCUSSION AND CONCLUSIONS 

The initial 7 and 8 TeV runs of the LHC have launched an impressive exploration of the 
vast parameter space of the MSSM and its extensions, yielding the most sensitive searches 
to date and in many cases already probing up to TeV mass scales. In this paper, we 
have sought to initiate a systematic extension of this program into the similarly vast realm 
of SUSY with non-prompt hnal-stage decays, as occur in a number of common and well- 
motivated scenarios within the MSSM such as R-parity violation, gauge mediation, and 
mini-split spectra. This has been accomplished by recasting seven existing searches for 
stable charged particles and highly displaced decays, and combining these with prompt 
searches. Our present survey has focused mainly on signals containing a sizable fraction of 
hadronic decays, including “natural” spectra with light stops and Higgsinos. The overarching 
conclusion is that, while very few long-lived particle searches are explicitly designed to be 
sensitive to these signals, the extremely low backgrounds and reasonably high acceptances 
of those searches nonetheless allow us to place tight limits. Indeed, we have typically found 
large patches of parameter space where multiple distinct search channels overlap. That said, 
we have identihed several places where signihcant improvements might still be possible. 

We hrst list here some conclusions of our investigations regarding the performance of the 
searches: 
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• In the long-lifetime limit, several searches have been carried out for stable charged 
R-hadrons. They are usually still sensitive down to cr ~ meter, catching rare late- 
decaying particles, and providing substantial overlap with explicit displaced decay 
searches. Limits on squarks and gluinos extend up to about 1 TeV. However, there is 
of course no sensitivity to long-lived neutral particles, such as (N)LSP Higgsinos. 

• The CMS displaced dijet search is extremely effective for essentially any decay topology 
involving any number of energetic quarks and gluons (including one or three, as well 
as many decays with leptons dl). It is almost universally the most powerful displaced 
decay search when hadronic decays dominate, including decays involving weak bosons. 
For CT ~ 10 cm, (stop) squark pairs are probed up to almost 1 TeV, gluinos typically up 
to above 1 TeV, and Higgsinos up to 600-800 GeV. Though the trigger px requirements 
are relatively harsh, good sensitivity is maintained for superparticles with strong or 
electroweak production cross sections down to 100 GeV mass, by picking up events on 
the high-pr tails. It seems likely that similar search strategies, perhaps capitalizing on 
different sets of triggers and displaced vertex criteria, could expand the model reach. 
In particular, within our own set of models the original dijet requirement reduces 
by 0(10-100) the acceptance for “monojet” decay topologies such as from q —)■ qG 
and g —)■ gG in GMSB. A standard jets-|-.^T style trigger could much more efficiently 
pick up this signal, and a single displaced-jet requirement would likely eliminate the 
remaining background. It would also be very useful to investigate this style of search for 
more traditionally motivated signals such as r —)• tG, which would require accepting 
low-multiplicity/low-mass vertices. For this signature and many others, it is important 
to understand what happens when the displaced particle is charged and can leave its 
own track segment. 

• Existing ATLAS searches for displaced decays inside the calorimeters and muon cham¬ 
bers should in principle supplement the region cr ~ few meters, where both the 
stable charged track and displaced dijet searches are becoming weaker. However, 
these searches were highly optimized for low-mass pseudoscalar signals, and achieve 
relatively limited sensitivity for displaced SUSY models. It seems likely that these 
searches could be modified to better capitalize on the energetic signatures of super¬ 
particles with lOO’s of GeV mass, where they are anyway most needed to extend the 
global search reach beyond cr ~ 1 m, especially for long-lived neutral particles where 
stable charged track searches are unusable. This could possibly be done using existing 
data from the specialized ATLAS HGAL and muon triggers. Such modified searches 
would still need to achieve very high efficiency in order to become competitive with 
the other searches, possibly benefiting from a single-candidate mode rather than their 
standard double-candidate. 
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• We have also considered three searches that capitalize on relatively clean leptonic 
signatnres: a displaced mnon in association with a displaced tracker vertex (not nec¬ 
essarily geometrically overlapping), a displaced dilepton pair, and generic e/r -1- X sig¬ 
natnres where the electron and mnon are each displaced. Interestingly, the /i-|-tracks 
search shows nontrivial sensitivity even to fnlly hadronic decays, provided that they 
contain energetic bottom qnarks. However, within the scope of models stndied here, 
all of these searches trnly become relevant for signals involving weak bosons, snch as 
stop and Higgsino decays in GMSB. The displaced dilepton search in particnlar offers 
improved sensitivity relative to displaced dijets when Z bosons are available, since the 
former can be constrncted with both lower impact parameter thresholds and lower 
Pt thresholds. The /x-|-tracks is also highly competitive, thongh somewhat weaker in 
lifetime coverage dne to geometric restrictions on the vertexing.^® 

• A variety of prompt searches have been directly applied to the decay topologies that 
we consider here. It is, however, mostly nnclear how effective these searches are when 
the decays become appreciably displaced, especially for searches involving leptons or 
traditional b- and r-jets. A first analysis in this direction was performed by ATLAS [10] 
for glninos in mini-split SUSY, indicating an approximately logarithmic degradation 
of mass reach with increasing lifetime (presnmably stemming from the onset of a 
linear falloff in displaced decay acceptance when cr > 0(m)). For searches at CMS, 
which rely mnch more on tracking and vertexing in the constrnction and validation 
of particle-flow jets, more significant degradations might be expected. While we have 
shown that CMS’s dedicated displaced dijets search can mnch more efficiently pick np 
the glnino signal in the lifetime range where ATLAS shows resnlts, we expect that there 
wonld be additional benefits to exploring other searches involving “many 6-jets” and/or 
“many r-jets” (possibly pins ^t), ideally with some allowances for nncharacteristically 
high displacements. Snch searches wonld (or perhaps already do) bridge the possible 
weakening in fnlly hadronic coverage aronnd 0(mm) lifetime stemming from CMS’s 
500 /im displaced dijet impact parameter cnt. As noted, the benefit of searching 
down to smaller nonzero displacements was already made clear in onr CMSB Higgsino 
results, where displaced dileptons was able to push the lifetime reach down by as much 
as an order of magnitude. 

We next discuss some of the physics implications of our findings: 

• Natural supersymmetry with light, promptly-decaying stops has been coming under 
progressively more pressure from a series of LHC searches. Moving to scenarios with 


It could be quite interesting to recast this search (as well as the displaced dijets) for the i — >■ leptonic 

RPV model, which was searched for in the e^ + X channel in [T^ . There are likely other models we have 
not touched upon for which this search could be uniquely sensitive. 
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displaced stop decays into jj via baryonic RPV or in GMSB, we are apparently 
forced into to even more nnnatnral regions of model space. This is particularly true 
in baryonic RPV, where prompt decay limits are currently very weak. However, very 
small A" and displaced stop decays are actually favored by cosmological arguments, as 
even X" > 10“® with sub-TeV stops would efficiently wipe out the baryon asymmetry 
of the universe in a standard thermal history with baryogenesis no lower than the 
weak scale [HH] • R is therefore becoming very difficult to simultaneously accommodate 
naturalness, baryogenesis, and LHC direct searches in such a scenario. In GMSB, 
where an NLSP stop is rather nongeneric but has been an interesting logical possibility 
for some time [83], the region with < rrtt is almost guaranteed to yield displaced 
decays (provided a/F > 10 TeV). The mi < mt possibility is now fully closed, as is 
essentially all model space below 500 GeV. 

• Natural supersymmetry with light Higgsinos is traditionally very challenging to probe 
via direct electroweak Higgsino production, though some limits are becoming available 
in general GMSB models where decays into Z bosons are appreciable [93] . In the 
presence of baryonic RPV, the multijet decays of the lightest Higgsino yield a very 
striking and highly constrained displaced signature, and furnish the only extant LHC 
direct production limits in that topology. The cosmological implications are more 
model-dependent, but again tend to disfavor the as-yet unprobed prompt decays [HH]- 
RPV Higgsinos with cr > 10 m would effectively act stable, and again become very 
difficult unless, as usual for neutralinos, they are produced in the decays of heavier 
colored superparticles and appear as ^t- In the GMSB case, as noted above, almost 
all of the searches that we have recast become sensitive, the only exception being 
stable charged particles. The mass/lifetime coverage is qualitatively similar to the 
RPV case, though the smaller fraction of visible energy, the typically smaller number 
of hard partons, and the smaller branching fractions into individual hnal states all 
contribute to slightly lower mass reach. Assuming that a natural Higgsino mass must 
be roughly below 400 GeV (naively corresponding to less than 5% hne-tuning of the 
Higgs boson mass), the displaced GMSB decay searches can probe a large fraction of 
the available space with a/F between 10 TeV and a several 1000 TeV. In particular, 
for ~ few 100 GeV, there is now a fairly hrm constraint a/F > 1000 TeV, unless 
the decay is dominated by hG. 

• Mini-split SUSY with sub-TeV gluinos is close to being fully ruled-out for any squark 

mass scale. An immediate escape hatch is to compress the spectrum to < 

100 GeV. We have also not studied in detail the limits on decays involving top quarks, 
though many other searches then open up, and it would be surprising if the limits 
become appreciably weaker. Non-minimal decay topologies involving electroweak-ino 
cascades could also be interesting to study, but would likely only yield signals that are 


35 


even more visible to displaced decay searches. 


• Generic colored snperparticles as (N)LSPs in either RPV or GMSB are to large degree 
rnled ont for any lifetime if the mass is below abont 1 TeV, again with the glnino limits 
tending to be several 100 GeV stronger than squark limits. Light RPV sqnarks with 
prompt decays would also face direct search difficulties similar to stops, but partially 
compensated by the higher multiplicity of flavor/chirality states. 


What else remains to be done? Within the context of RPV (both baryonic and leptonic), 
a more thorough survey of the current status of different LSPs and flavor structures along the 
lines of |99] and [28l 1100] seems warranted. Leptonic RPV in particular has a quite large set 
of possible couplings. Spectra with “electroweak” LSPs besides Higgsinos, namely sleptons 
or gauginos, also deserve further attention, as they can become much more visible than they 
would be if their decays were prompt. For general gauge mediation, we have emphasized in 
Section III B that the full set of possible NLSPs is (rather remarkably) almost fully covered. 
The major exceptions are again sleptons and winos, with the latter offering an interestingly 
varied array of different signatures. We again expect that the existing set of displaced and 
prompt searches have much to say about all of the above models, though in many cases 
coverage may still be entirely lacking, unnecessarily weak, or ambiguous given the current 
limitations of making public the general analysis acceptances. 

With the upcoming Run 2 of the LHG, the mass reach for the models that we have 
explicitly studied might be expected to roughly double, assuming that similar analyses will 
be undertaken. We encourage the experiments to continue their displaced decay search 
programs at an even greater level of breadth so that interesting signals are not left behind. 
We also hope that future recasts are better facilitated by more explicit discussions of analysis 
acceptances, less tied to one or two specific fully simulated models in limited kinematic 
ranges. Endeavors like ours should ideally not require as much from-scratch calibration, 
extrapolation, and guesswork, as detailed here in Section [IT] and in the Appendices below. 
Of course, the need to facilitate more general model interpretations becomes even more 
pressing if a discovery is made. Works along the lines of GMS’s stable charged particle 
efficiency maps jH] are a step in the right direction. But we emphasize that even coarse 
parametrizations such as the ones that we have developed can prove invaluable, especially 
if directly compared against full internal simulations by the collaborations. We hope that 
our work, which has further clarified the extreme power and broad model reach of these 
searches, spurs further activities in these directions, and we look forward to the next round 
of LHG displaced search results. 
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Appendix: Detector Simulations and Calibrations 

The set of recastings performed in this paper are all based on a small set of custom 
detector simulation codes. These are meant to capture the main features relevant to each 
analysis while bypassing the many highly complex details of the real detector response. 
Since the displaced SUSY models studied here have never appeared in public collaboration 
results, calibration of our simulations must rely on the specihc physics models that appear 
in the experimental literature. Our aim has been to build simulations that reproduce the 
known experimental acceptances for new physics at the 0(1) level, which we are often able 
to accomplish even with very minimal treatment of detector issues. In most cases this level 
of agreement is adequate to define reasonable estimates of the true sensitivity contours in the 
mass-lifetime plane, as the sensitivity can be an extremely steep function in both variables. 
However, typically we can achieve even 0(10%) agreement with the experimental results, 
either “out of the box” or by adjusting ad hoc efficiency factors. A single tuning of the latter 
type is often able to reproduce the results for a broad range of models. 

Below, we provide complete descriptions of our detector simulations and their calibrations. 


A. CMS Heavy Stable Charged Particles 

The simulation required for the CMS heavy stable charged particle search [Sj can be 
quite minimal, since CMS searches for the same physics signal that we do: long-lived stop 
and gluino R-hadrons. The only novelty that we introduce is the finite lifetime, which 
effectively adds an additional factor to the overall R-hadron acceptance. Decays must occur 
completely outside of the muon system, with no visible particles from the decay pointing 
back. We define the outer edge of this active volume simply as a cylinder with 8 m radius 
and 10 m half-length.^® 


Reference [9] recommends using a slightly smaller radius of 7 m and a slightly longer half-length of 11 m. 
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Calibration is trivial, as we define onr analysis in Section IIA to reproduce the results 
of [8] for long lifetimes. 


B. CMS Displaced Dijets 

The CMS displaced dijet search [miss] utilizes that detector’s highly precise tracking 
capabilities. Nonetheless, the tracker is not perfect, especially in the extreme cases of 0(m) 
displacements. At such large displacements, particles may not traverse enough detector 
layers to furnish a reconstructable string of hits, or may lead to hit patterns that not usable 
given a detector geometry and track-finding software that is highly optimized for tracks 
originating very close to the beampipe. In particular, CMS observes a dramatic drop in 
tracking efficiency versus vertex radius in simulation |ini] , falling to zero beyond 60 cm. The 
tracker also becomes highly inefficient for soft particles, especially for transverse momenta 
that are small enough for the particles to spiral-out. The efficiency experiences a rapid 
turn-on near 1 GeV. 

It is not possible to reproduce the tracking and vertex-finding performance in detail 
here. Instead, we rely on simplistic parametrizations, and validate them against the results 
presented in the analysis note |56]. As a “zeroth-order” approximation, we could consider 
the tracker to be perfect within some fiducial volume. Our most naive version applies hard 
cutoffs at r = 60 cm and pt = 1 GeV. Based on the discussions in [IHIEE], and Fig. 3 in [T8] . 
we also include a corresponding longitudinal position cutoff at z = 55 cm, as well as a cutoff 
in transverse impact parameter at 30 cm. However, we have found by comparing to [HS] 
that this zeroth-order treatment is far too idealized. GMS’s new physics models with 0(cm) 
lifetimes are reproduced fairly well, with reconstruction rate estimates already typically 
within about 25% of GMS’s. (This result by itself indicates the amazing performance of 
the GMS tracker system for charged particles originating near the beampipe.) But the 
efficiencies for the longer 0(m) lifetimes come out too large by up to a factor of 3.5, and 
distributions of variables such as the reconstructed vertex radius often significantly disagree. 

We therefore use a slightly more aggressive parametrization for our nominal detector 
simulation. We continue to apply a hard cutoff at = 1 GeV. For the geometric limitations, 
instead of hard cutoffs, we apply a nontrivial track-finding probability. We construct this 
probability by starting off with a prompt tracking efficiency of 90%, and multiplying it by 
the product of three linearly-falling probability factors, one each for r, z, and transverse 
impact parameter. Each factor starts at unity, and falls to zero at the geometric cutoff. 
Finally, in addition to these track-by-track efficiencies, we apply an efficiency for identifying 
each candidate vertex. This efficiency is the square of a linear falloff function versus radius, 
intercepting zero at 60 cm. The combination of these ad hoc efficiency factors significantly 
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CMS Displaced Dijets 




FIG. 11; A comparison of our detector simulation to the CMS displaced dijet analysis, illustrating 
the ratio of individual displaced decay reconstruction rates to CMS for light quark decays (left) 
and the ratio of ratios between 6-quarks and light quarks (right). The 14 benchmark models 
are labeled by scalar / pseudoscalar masses in GeV and by cr in cm. On the left plot, we show 
our nominal simulation (black), our nominal simulation without the vertexing efficiency factor 
(pink), and our “zeroth-order” simulation with perfect tracking within the fiducial volume (red). 
Error bars are monte carlo statistics from our simulations. Grey bands indicate CMS’s efficiency 
uncertainties. 

improves our predicted rates and distributions relative to CMS.^^ 

For the rest of the detector, we use a very minimal treatment to capture the essential 
elements of the geometry. The tracker is surrounded by imaginary surfaces corresponding 
the ECAL and HCAL faces. The ECAL face corresponds to a closed cylinder of radius 1.3 m 
and half-length 3 m. The HCAL face corresponds to a closed cylinder of radius 1.8 m and 
half-length 3.75 m. Charged particles that fail tracking due to the above inefficiencies are 
analytically propagated through the 3.8 T magnetic held to the appropriate calorimeter face: 
ECAL for electrons, HCAL for charged hadrons (failed muon tracks are discarded). Photons 
and neutral hadrons are propagated to the appropriate calorimeter face along straight lines. 
Particles are absorbed at the calorimeter face, and replaced by anonymous massless pseu¬ 
doparticles carrying the original particle energy, but with momentum oriented to toward the 


The additional ad hoc vertexing efficiency causes a more rapid degradation of reconstruction rates versus 
radius than what is reported by CMS in their supplementary online materials |102j . However, an alterna¬ 
tive approach of convolving the reported rates versus radius with our simulated number of decays versus 
radius similarly disagrees with the CMS analysis note [56]. We do not attempt to resolve this apparent 
discrepancy, but haven chosen to parametrize our simulation to reproduce the integrated rates of the note, 
rather than the differential rates of the online material. This is in any case the more conservative choice. 


39 



















impact point.For successfully reconstructed tracks, we use the momentum vector at their 
extrapolated point of closest transverse approach to the beamline (where the 2D impact 
parameter is dehned). Track and calorimeter energy/momentum measurements are treated 
as perfect, since the effects of energy smearing are subdominant to our other modeling un¬ 
certainties. The tracks and calorimeter pseudoparticles are the inputs into our jet clustering, 
excluding any identihed isolated electrons or muons as dehned in the next subsection. 

We plot our calibration results of the per-decay reconstruction rate (“acceptance times 
efficiency” in CMS’s language) relative to CMS in Fig. 11 These comparisons are based 
on the set of Hidden Valley [103] inspired simplihed models used by CMS, consisting of a 
heavy Higgs-like scalar, produced through gluon fusion, that decays into a pair of long-lived 
pseudoscalars. The pseudoscalars then undergo displaced decays into dijets. (We have not 
undertaken calibrations against the —)• jjijj signatures that are studied in the more recent 
preprint Cl-) The steep slopes for our zeroth-order detector’s reconstruction rates rela¬ 
tive to CMS indicates that simulation’s failure to correctly account for inefficiencies at high 
displacements. A remnant of this slope remains after applying our tracking efficiencies but 
before applying our vertexing efficiencies. The nominal detector is generally in agreement 
with CMS for all lifetimes to within 20%, with the notable exception of the 200 GeV scalar 
model. We do not consider this a serious issue, as this model has high sensitivity to initial- 
state radiation modeling. The other models exhibit a good spread of overall masses, mass 
hierarchies, and lifetimes. In particular, the 400/50 and 1000/150 models tend to produce 
dijets near the edge of the jet clustering radius, and we see that we tend to slightly underes¬ 
timate their reconstruction rates relative to models with widely-separated jets. Nonetheless, 
the highest dijet mass covered is 350 GeV, whereas some of our SUSY models go above 
1 TeV. We assume that there are no dramatic changes in efficiencies as we scale up in 
mass, though a broader set of simulated models from CMS would help to clarify the actual 
behavior. 

Fig [TT] also indicates our ability to reproduce reconstructions with heavy flavor, by show¬ 
ing a double-ratio of reconstruction rates. The numerator is our estimated ratio of rates for 
A —)■ 66 relative to X decays to light flavors. The denominator is CMS’s estimate of the 
same ratio. For this analysis, we pretend that the secondary displacements from the bottoms 
either cannot be resolved or are effectively ignored by the adaptive vertex hnder. The agree¬ 
ment is generally seen to be quite reasonable, with a handful of outliers disagreeing at more 
than 20%. The fact that we achieve such good agreement without separately displacing the 


For photons and electrons produced from displaced vertices in between the ECAL and HCAL faces, the 
energy is deposited without propagation, and the momentum vector pointed towards the vertex. Any 
particle produced within the HCAL body is also absorbed in this way. For this purpose we assign an 
outer HCAL surface of radius 2.8 m and half-length 5.5 m. Particles produced outside of this volume are 
ignored. 
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bottoms strongly suggests that the dominant differences in light flavor and heavy flavor effi¬ 
ciencies stems from the different visible particle multiplicities and kinematics. There remains 
a question of whether small primary displacements at the mm-level may have resolvable sec¬ 
ondary displacements according to CMS’s adaptive vertex fitter, such that only the “best” 
of the truth vertices from each decay actually contributes. (E.g., i ^ bb could produce up 
to three separate vertices, one from the radiation before the 6’s hadronize, and two more 
from the 6-hadron decays.) Presumably, this would degrade the efficiency, since fewer tracks 
would be usable from any given decay. However, CMS does not provide enough information 
to infer exactly what happens for primary displacements below the cm-scale. In the main 
results in our paper, we only consider the extremely conservative assumption that vertices 
with bottom/charm quarks and sub-cm displacements experience complete reconstruction 
failure. 

Finally, we comment on our modeling of the variables that go into CMS’s multivariate 
vertex/cluster discriminant variable. The discriminant is defined as a ratio of products of 
normalized p.d.f.’s over four variables: vertex track multiplicity, vertex positive IP fraction, 
cluster track multiplicity, and cluster RMS. The “cluster” is formed as described in Sec¬ 
tion IIB We plot our predicted Hidden Valley signal distributions for the four variables and 
the multivariate discriminant in Fig. 1^, along with the CMS predictions. For both our own 
simulation and CMS’s, we form the discriminant using CMS’s p.d.f.’s. (The signal discrim¬ 
inant distribution that would be predicted by CMS is generated by us through toy monte 
carlo from the individual variables’ histograms.) The agreements are mostly reasonable, 
though there can be differences in the tails. The fraction of events passing the discriminant 
cuts are nonetheless nearly identical to CMS. The signal’s discriminant distributions at high 
values appear to be mostly driven by the multiplicity variables, which we model relatively 
well. 


C. CMS Displaced Dileptons 

For the CMS displaced dilepton search [IS], we use the same detector simulation discussed 
in the previous subsection, though without the vertex reconstruction penalty as a function 
of radius. We also apply the pr-dependent lepton ID efficiencies provided in the appendix 
of |104j . though we divide these by 0.9 to approximately deconvolve the prompt track-finding 
efficiency which we have already accounted for. 

Similar to the displaced dijet search, the benchmark models feature a Higgs-like scalar 
decaying to a pair of long-lived pseudoscalars, but with the latter decaying to or 

■ Unlike the displaced dijet search, detailed tables of acceptances and efficiencies are 
not provided. However, a handful of specific numbers are given, and the observed cross 
section limits can also be used to infer overall reconstruction rates. 

For the 1000 GeV / 150 GeV models, GMS’s individual pseudoscalar reconstruction 
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FIG. 12: A comparison of our detector simulation to the CMS displaced dijet analysis, illustrating 
expected reconstructed dijet counts versus several discriminator variables (Fig. 1 of [56]). Con¬ 
tinuous histograms indicate CMS predictions for QCD background (blue) and different scalar / 
pseudoscalar masses: 1000 GeV / 350 GeV with cr = 35 cm (red), 400 GeV / 150 GeV with 
CT = 40 cm (green), and 200 GeV / 50 GeV with cr = 20 cm (black); all with cr —)• 10 /rb for 
visualization. Data points with error bars are our simulation predictions, and are color-matched 
and normalized to the corresponding CMS model histogram. The multivariate discriminant cuts 
for High-Lxy and Low-Lxy signal regions are also indicated. 
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CMS Displaced Dilepton: H(1000)-»4e 


CMS Displaced Dilepton: H(200) ^ 4p 




FIG. 13: A comparison of our detector simulation to the CMS displaced dilepton analysis by 
inverting their reported limits, illustrating the approximate ratio of event reconstruction rates to 
CMS for 1000 GeV scalar cascade decays to electrons (left) and 200 GeV scalar cascade decays to 
muons (right) versus pseudoscalar lifetime (based on the limits presented in Figs. 4 and 5 of |18jl. 
On the left plot, we show pseudoscalar masses of 20 GeV (pink), 50 GeV (blue), 150 GeV (red), 
and 350 GeV (black). On the right plot, we show pseudoscalar masses of 20 GeV (red) and 50 GeV 
(black). Error bars are monte carlo statistics from our simulations. All models are evaluated at cr 
in powers of 10, but are offset slightly horizontally for clarity. 

rates for cr = 1 cm in the electron (muon) decay channel are given as 36% (46%). Our 
predictions are in decent agreement, at 36% (40%). For cr = 20 cm in the electron (muon) 
decay channel, CMS gives 14% (20%). Our predictions here are 9% (10%), indicating too-low 
reconstruction rates by 0(1). This is likely due to our ad hoc track-hnding efficiency’s linear 
falloff being too steep for this analysis. Indeed the supplementary online material of |105j 
indicates fairly stable reconstruction efficiency for this model out to r ~ 50 cm, where our 
simulation would predict nearly zero. It seems quite likely that the low track multiplicity of 
the decay contributes to a higher rate of successful displaced track reconstructions, relative 
to that of the displaced dijets. However, we conservatively continue to use our tracking 
efficiency factors derived for the latter analysis, especially as in some SUSY models the 
lepton pair may be produced in association with hadronic tracks from the same vertex. This 
does not lead to a runaway loss of efficiency at high lifetime for simple dilepton decays, as 
we will see. 

The calculation of limits from the signal reconstruction rates is in principle nontrivial, 
involving incorporation of backgrounds and various sources of systematic errors. However, 
given that the analysis is ultimately both background-free (as inferred from a control re¬ 
gion) and has zero observed events, and that the naive Poissonian 95% C.L. limit for such 
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CMS Displaced Electron and Muon 



FIG. 14: A comparison of our detector simulation to the CMS displaced electron and muon analysis, 
illustrating the ratio of individual displaced decay reconstruction rates to CMS for the 500 GeV stop 
benchmark model. The bins are labeled by analysis signal region (SRI for lower impact parameters, 
through SR3 for higher impact parameters) and by the stop’s cr in cm. Grey bands indicate CMS’s 
efficiency uncertainties, with the stop cross section uncertainty quadrature-subtracted. 


an experiment is ~ 3 signal events, we can perform a back-of-the-envelope estimate of the 
reconstrnction rates. The explicit CMS per-candidate rates discnssed above are reprodnced 
at the 10% level nsing this method, providing a good cross-check. Fig. 1^ illustrates our es¬ 
timated relative rates for a few example models. The general behavior is that our simulation 
appears to be less efficient than reality, especially at very low and very high lifetimes. The 
former is perhaps unsurprising given our coarse modeling of the impact parameter cut, and 
the latter could be due to our overzealous tracking efficiency falloff. Nonetheless, some of 
the model points still appear to exceed the CMS rate, even at higher lifetimes, motivating us 
to keep these somewhat artihcial effects to help prevent us from inferring too-strong limits 
on SUSY models. 


D. CMS Displaced Electron and Mnon 

For the CMS displaced electron and muon search [15], we continue to use the detector 
simulation described in Appendix To better match the efficiencies reported by CMS in 
the present analysis, we apply a flat event-by-event weight of 0.80. 

The physics model studied by CMS is RPV stop pair production, but with each stop 
decaying into bl~^, with equal branching fractions into each of the three lepton generations. 
It is not explicitly stated in this analysis whether the displaced leptons could suffer any loss 
of efficiency if produced from the decay of a charged stop-hadron, leading to an “exploding 
track” topology rather a displaced vertex with no string of hits tracing back to the primary 
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vertex. In any case, since the analysis focuses on impact parameters below 2 cm, it is 
mostly sensitive to decays that occur before reaching the pixels. We treat charged and 
neutral stop-hadrons identically, assuming that this is not a major issue. CMS gives explicit 
reconstruction rates for a 500 GeV stop at lifetimes of 0.1 cm, 1 cm, and 10 cm. Our detector 
simulation does a good job of reproducing all of these numbers to 20% accuracy, as indicated 


in Fig. 14 


E. ATLAS Muon Spectrometer 


The ATLAS muon spectrometer search [20] is again very difficult to model without access 
to both a full ATLAS detector simulation and the exact reconstruction algorithms. Here, 
we simply parametrize all of these with hxed efficiency factors. Our simulation dehnes an 
active trigger volume within r = [4.0, 6.5] m and |? 7 | < 1.0. Displaced particles that decay 
in this region are given a 50% chance to hre the muon Rol cluster trigger. For events 
passing the trigger, two displaced decays must be reconstructed in either the muon barrel 
or muon endcap, respectively dehned as the volumes r = [4.0, 7.5] m and \r]\ < 1.0, or 
\ri\ = [1.0, 2.5] and z = [8,14] m. Both of these decays must also occur within the data 
acquisition timing window, designed to follow particles moving through the detector at the 
speed of light. We choose a maximum delay of 7 ns, which we have independently calibrated 
to the extremely timing-sensitive ATLAS 120 GeV/ 40 GeV model results. Our calibrated 
choice of maximum delay indeed corresponds to the end of the efficiency plateau versus time 
for the trigger |106j . We assign each displaced vertex that survives this cut a reconstruction 
rate of 40%, irrespective of which ones were capable of bring the Rol trigger. Finally, in 
order to very approximately account for possible isolation failures when a displaced decay 
points back to the detector volume, we limit the amount of visible transverse energy bowing 
back into the HGAL to 15 GeV. For this purpose, we debne the outer surface of the HGAL 
as a cylinder of radius 4.25 m and half-length 6 m. This cut tends to have only modest 
ebect on ATLAS scalar models, but could become important for our SUSY models at higher 
mass. 

Similar to the GMS displaced dilepton analysis above, we calibrate against the ATLAS 
models by extracting their overall efficiencies from the reported limits. Given that this 
analysis is background-free, the limits are assumed to correspond to ~ 3 signal events. The 


result of this comparison is shown in Fig. lA We are able to reproduce all of the ATLAS 
results to within about 20%. (The one outlier point occurs on a very steep part of the 
efficiency curve.) 
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FIG. 15: A comparison of our detector simulation to the ATLAS muon chamber analysis using their 
reported limits, illustrating the ratio of event reconstruction rates to ATLAS versus pseudoscalar 
lifetime (based on the limits presented in Fig. 3 of [20]). We show scalar / pseudoscalar masses of 
120 GeV / 20 GeV (solid black), 120 GeV / 40 GeV (dashed black), 140 GeV / 20 GeV (solid red), 
and 140 GeV / 40 GeV (dashed red). Error bars are monte carlo statistics from our simulations. 
All models are evaluated at common cr of (0.5,1,1.5,2,3,4,5,10,15,20) m, but are offset slightly 
horizontally for clarity. 

F. ATLAS Low-EM Jets 

For our SUSY models to pass the ATLAS low-EM jets search ng, both displaced 
decays must occur within an effective HCAL volume: a “barrel” with \r]\ < 1.7 and 
r = [2.25,3.35] m; or an “endcap” with \ri\ < 2.5, r < 2.0 m, and z = [4.25,5.0] m. 
To approximate ATLAS’s isolation cuts against activity in the ECAL and tracker, we ap¬ 
ply a flat efficiency factor of 0.55 per displaced decay.(As usual, events with one or two 
charged displaced particles, which would leave high-p^ tracks, are not considered.) 

Our detector simulation assumes perfect and immediate absorption of the visible decay 
energy within our active HCAL volume. However, we conservatively veto events where any 
visible hnal-state particle from either displaced decay points back towards the ECAL volume, 
which we take to be r < 2.0 m and \z\ < 4.0 m. The effect of this veto is modest for most of 
atlas’s Hidden Valley models, but notably has a nearly 50% effect on the acceptance for 
the highest mass models, actually improving our agreement (see below). The cut is especially 
relevant for our SUSY models, which generally contain many more decay particles produced 
at large angles. Practically, for us the cut serves as a proxy for any number of unknown 

ATLAS’S reported efficiencies are smaller than this. However, those also fold in the efficiencies for passing 
the pt cuts within their studied physics models. 
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FIG. 16: A comparison of our detector simulation to the ATLAS low-EM jet analysis, illustrat¬ 
ing expected event counts versus pseudoscalar lifetime (Fig. 4 of [IH]). Curves indicate ATLAS 
predictions for different scalar / pseudoscalar masses: 126 GeV / 10 GeV (red solid), 126 GeV / 
25 GeV (black dashed), and 140 GeV / 40 GeV (pink dotted). Data points with error bars are our 
simulation predictions, and are color-matched to the corresponding ATLAS model curve. 

requirements on the energy pattern and timing in the HCAL, in addition to the explicit 
limits on the nearby ECAL activity. As such, it serves as our largest source of modeling 
uncertainty on this analysis. 

To validate our detector simulation and calibrate the efficiency factor, we compare to AT¬ 
LAS’S physics models, which again feature singly-produced Higgs-like scalars that promptly 
decay into a pair of pseudoscalars with displaced decays. The pseudoscalars dominantly 
decay into heavy quarks, with an 0(10%) branching ratio to r-leptons, though we use a 
simpler model with decays only to light quarks. (Our model therefore does not include 
energy losses to neutrinos, which we do not expect to be a major effect.) We have found 
good agreement with the reconstruction efficiencies of the eight benchmark models, detailed 
in the appendix of [I9], though for 100 GeV scalars our efficiencies come out 0(1) smaller. 
In any case, these models barely pass the Et cuts, and are signihcantly more sensitive to 
the detailed turn-on of the efficiency with scalar mass. We have also studied the efficiency 
versus pseudoscalar lifetime for three of the mass points. The agreement is illustrated in 
Fig. [T^ and is typically at the 10% level. 


G. ATLAS Displaced Muon Plus Tracks 

Similar to our simulation of the CMS displaced dijet analysis, for the ATLAS muon plus 
tracks search HI we assign each track a reconstruction probability. We ignore all tracks 
from vertices with r > 18 cm or ^ > 30 cm, as well as cylindrical regions near the beampipe 
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FIG. 17: A comparison of our detector simulation to the ATLAS muon plus tracks analysis, 
illustrating expected event reconstruction efficiencies versus neutralino lifetime (Fig. 4 of |14j). 
Colored bands indicate ATLAS predictions for different squark / neutralino masses: 700 GeV / 
494 GeV (blue), 700 GeV / 108 GeV (red), and 1000 GeV / 108 GeV (purple). Data points 
with error bars are our simulation predictions, and are color-matched to the corresponding ATLAS 
model curve. 


or the pixel layers: r = [2.5, 3.8] cm, [4.5, 6.0] cm, [8.5, 9.5] cm, and [12,13] cm. Within 
the active tracking volnme, we nse a tracking efficiency of 0.85 x (1 — r/(24 cm)) x (1 — 
^/(30 cm)). For mnon identihcation, which reqnires a very high-qnality inner track, we 
apply an additional hxed efficiency factor of 70%. We apply another hxed factor of 70% 
for snccessfnlly matching a mnon to a displaced vertex. These ad hoc factors reprodnce 
the qnalitative vertex reconstrnction efficiency behavior versns (r, z) in Fig. 3 of na , as 
well as the ratio between event selection efficiencies before/after the mnon-vertex matching 
reqnirement in Table 5 of that note.^° The latter is particnlarly relevant for the version of 
the analysis that we rnn for onr recasts, which does not reqnire the mnon-vertex matching 
nsed in the nominal analysis. 

ATLAS stndies a set of SUSY models with leptonic RPV, where sqnark pairs promptly 
decay into jets pins long-lived LSP nentralinos, which then each nndergo a displace decays 
into a mnon pins two qnarks. Three mass points are considered: MH (700 GeV / 494 GeV), 
ML (700 GeV / 108 GeV), and HL (1000 GeV / 108 GeV). Fig. 17 shows the event recon¬ 
struction efficiencies versus lifetime for these models after full selection cuts, as predicted 
by our simulation and by ATLAS’s full simulation. (The individual vertex reconstruction 
efficiencies are approximately half of the event reconstruction efficiencies, and show nearly 


The table suggests that the vertex reconstruction efficiencies, before track multiplicity and vertex mass 
cuts, are close to one. However, matching between outer muons and inner tracks, as well as between 
muons and vertices, exhibits these nontrivial efficiencies. 
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identical behavior.) The large gap in efficiencies between MH and ML/HL is reproduced, 
as is the shift in the peak versus lifetime and the falloffs at low and high lifetimes. How¬ 
ever, the ~ 30% gap between ML and HL is not reproduced, except at higher lifetimes, 
where the higher boost of HL tends to feel the tracking efficiency losses at large radii. (The 
higher boost would also lead to smaller impact parameters at a given decay radius, but this 
is largely offset by the time-dilated decay length.) The difference in modeling is possibly 
attributable to the fact that our simulation does not account for how displaced tracking and 
vertexing efficiencies change with track density, nor to possible issues in impact parameter 
reconstruction at small angles, all of which the lighter neutralino could be particularly sen¬ 
sitive to. (To get a sense of this sensitivity, deleting one track from each vertex would cause 
the efficiency to fall by 20% due to failures of the cut on the number of tracks.) There could 
also be effects on the global muon reconstruction. Nonetheless, the size of the mismodeling 
is below 0(1), and appears to be mainly relevant for relatively low-mass decays at relatively 
high boost. 
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